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AND  TIME  AND  THE  INFLUENCE  OF  HABITAT 
COMPLEXITY  ON  THEIR  INTERACTIONS 


avoided  these  sites.  Hydrologic 


disturbance  resulted  in  a system  that  is  both  spatially  and 
temporally  dynamic  resulting  in  variable  assemblage 
structure  and  composition  at  affected  sites.  Distribution 
patterns  of  tadpoles  within  wetlands  were  related  weakly  to 
habitat  features.  Both  tadpole  and  aquatic  insect  predator 
numbers  varied  among  wetlands  and  across  time,  but  were 
distributed  relatively  evenly  within  wetlands.  Thus, 
selection  for  microhabitat  at  the  within-wetland  scale  did 
not  appear  to  be  strong.  Distribution  patterns  of  tadpoles 
appear  to  be  more  strongly  related  to  macrohabitat  selection 
of  breeding  sites  by  adults  than  by  microhabitat  selection 
within  wetland  sites.  Experiments  demonstrated  that 
increased  cover  can  reduce  predation  rates  by  aquatic 
insects  on  tadpoles  that  reduce  activity  in  response  to 
predators  (B.  utricularia) , and  those  that  do  not  IHyla 
squirella) . Thus,  although  distribution  patterns  in  the 
field  were  not  correlated  strongly  to  habitat  features, 
tadpole  success  still  may  be  influenced  strongly  by  them. 
Finally,  reductions  in  activity  in  response  to  predators  had 
growth  and  developmental  costs  in  R.  utricularia  that  were 
dependent  on  background  resource  level.  Tadpoles  on  low  food 
metamorphosed  at  later  dates  but  at  larger  sizes  when 
exposed  to  predators.  Tadpoles  on  high  food  also  had 
prolonged  development  but  were  smaller  compared  to  tadpoles 
not  exposed.  Thus,  assessment  of  responses  to  both  resource 


level  and  threat  of  predation  may  be  necessary 
understanding  the  consequences  of  cost-benefit 
mediated  by  behavior. 
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CHAPTER 


INTRODUCTION 

Overview 

One  of  the  major  goals  in  ecology  is  to  determine  the 
factors  that  influence  patterns  of  distribution  and 
abundance  of  species.  This  task  is  not  a simple  one  because 
many  factors  impinge  upon  the  ecological  and  evolutionary 
processes  that  result  in  observed  relationships  among 
species.  Resource  availability  and  distribution  change 
along  many  spatial  and  temporal  scales,  and  the  biotic 
interactions  among  species  change  similarly.  The  resulting 
heterogeneous  environment  offers  choices  to  the  organisms 
living  within  it.  How  an  individual  responds  to  habitat 
heterogeneity  may  influence  food  availability  and  predation 
pressure,  and  thus  can  have  important  influence  on 
population  dynamics  and  ultimately  evolutionary  processes 
(Holt  1987). 

Examination  of  the  factors  influencing  distribution 
patterns  of  many  anurans  is  complicated  by  their  complex 
life  cycles,  which  include  a terrestrial  adult  stage  and  an 
aquatic  larval  stage  (Wilbur  1980).  Examinations  of  habitat 
use  and  resource  partitioning  of  anuran  assemblages  have 
focused  mostly  on  the  adult  stage  (Bowker  and  Bowker  1979; 
crump  1971,  1982;  Lizana  et  al.  1990;  Toft  1985).  Analyses 


of  larval  assemblages  usually  address  temporal  partitioning 
(Di::on  and  Heyer  1968j  Weist  1981)  or  partitioning  relative 
to  permanent  versus  temporary  pond  sites  (Gascon  1991; 
Woodward  1983).  Host  authors  come  to  similar  conclusions 
regarding  habitat  partitioning  of  larval  assemblages; 
resources  are  partitioned  along  spatial  (temporary  versus 
permanent  ponds)  and  temporal  (breeding  season)  axes . less 
research  has  been  done  on  larval  food  partitioning;  however, 
evidence  suggests  that  partitioning  along  the  food  axis  is 
not  as  strong  as  that  along  temporal  and  spatial  axes 
(Duellman  and  Trueb  1986;  Heyer  1974,  1976;  Seale  and 
Beckvar  1980) . 

These  studies  provide  important  information  regarding 
possible  selection  pressures  relative  to  habitat  use  and 
resource  partitioning  (Crump  1982;  Magnussan  and  Hero  1991). 
However,  they  do  not  provide  a complete  picture  of  the 
spatial  ecology  of  tadpoles  because  they  focus  on  one  scale. 
Distribution  patterns  of  organisms  can  be  examined  on  a 
coarse-grained,  macrohabitat  scale  or  a Cine-grained, 
microhabitat  scale.  Factors  that  can  influence  tadpole 
success,  such  as  competitor  or  predator  abundance  and 
composition,  vary  in  space  not  only  among  wetland  siLes 
(Gascon  1992;  Smith  1983;  Werner  and  McPeek  1994;  Woodward 
1983)  but  also  among  patches  within  wetlands  (Banks  and 
Beebee  1988;  Diehl  1988,  1992).  Therefore,  the  distribution 
patterns  of  tadpoles  at  both  the  macrohabitat  (breeding 


site)  and  microhabitat  (patch  site)  scale  can  have  a large 
influence  on  tadpole  growth  and  survival,  and  thereby 
assemblage  structure. 

Although  the  importance  of  considering  scale  has  been 
widely  discussed  (Wiens  1989,  1992),  few  studies  of  tadpole 
ecology  have  addressed  this  issue.  Processes  that  are 
important  population  regulatory  factors  at  one  scale  of 
investigation  may  fail  to  operate  at  other  scales  (Wiens 
1992).  For  example,  predation  has  been  identified  as  an 
important  process  influencing  anuran  breeding  site  selection 
(Magnusson  and  Hero  1991).  However,  how  predators  and 
tadpoles  interact  at  the  within-wetland  scale  is  poorly 
understood.  Because  ecological  patterns  are  not  independent 
of  scale,  examination  of  spatial  and  temporal  distribution 
patterns  at  several  scales  should  provide  a more  complete 
understanding  of  the  importance  of  a process,  such  as 
predation,  in  regulating  tadpole  populations  (Wiens  1989, 


Patterns  of  Distribution  Among  Wetlands 
Patterns  of  distribution  of  larval  anurans  among 
aquatic  sites  are  largely  determined  by  two  related  factors,- 
the  permanence  of  the  site  (i.e.,  hydrology)  and  the 
composition  of  predators  (Shelly  I995j  Werner  and  McPeeh 
199<l),  These  factors  are  strongly  interrelated  because 
predator  composition  changes  along  the  hydrologic  gradient. 


Fish  predators  dominate  in  permanent  sites;  however,  systems 
that  have  a cycle  of  filling  and  drying  prevent  predatory 
fish  populations  from  becoming  established.  A switch  from 
fish-dominated  to  invertebrate  or  salamander-dominated 
predator  systems  occurs  along  the  gradient  so  that  tadpoles 
developing  in  temporary  and  ephemeral  sites  are  subject  to  a 
different  suite  of  predators. 

Because  biotic  and  abiotic  conditions  change  along  a 
hydrologic  gradient,  tadpoles  developing  in  different 
wetlands  face  a different  set  of  challenges.  Traits  that 
increase  developmental  success  at  one  end  of  the  hydrologic 
gradient  may  not  be  suitable  at  the  other  end.  For  example, 
interspecific  variation  in  anti-predator  mechanisms  is 
related  to  differences  in  breeding  habitat  (Asevedo-Ramos  et 
al.  1992;  Kats  et  al.  1988;  Woodward  1983).  The  larvae  of 
several  species  that  breed  in  permanent  sites  are 
unpalatable  to  fish,  whereas  species  that  breed  in  ephemeral 
sites  without  fish  do  not  possess  this  anti-predator 
mechanism  (Formanowicz  and  Brodie  1982;  Kats  et  al.  1988). 

In  addition,  behavioral  adjustments,  such  as  reduced 
activity,  are  often  absent  in  species  that  breed  in 
ephemeral  sites  (Kats  el  al . 1988;  Woodward  1983). 

Decreased  foraging  activity  may  reduce  growth  or 
developmental  rates  (Shelly  1992,  1995;  Skelly  and  Werner 
1990).  At  ephemeral  sites,  mortality  threats  from  predators 
are  often  not  as  significant  as  the  threat  from  habitat 


drying  (Newman  1987;  Rowe  and  Dunson  1993,  19951 . Thus, 
behavioral  traits  for  reducing  predation  risk,  while 
beneficial  in  permanent  habitats,  can  lead  to  increased 
mortality  risk  at  ephemeral  sites. 

The  trade-off  between  rapid  development  and  anti- 
predator responses  appears  to  be  a major  factor  limiting  the 
distribution  of  anurans  to  certain  types  of  breeding  sites 
(Skelly  1995;  Werner  and  McPeefc  199*11  . Weak  anti-predator 
responses  that  may  increase  developmental  rates  in  ephemeral 
sites  would  usually  result  in  mortality  in  permanent  sites 
(Skelly  1994;  Woodward  1983).  As  a broad  generalisation, 
species  that  breed  in  permanent  sites  are  passively  excluded 
from  ephemeral  sites  by  long  larval  developmental  periods, 
whereas  species  that  breed  in  ephemeral  sites  are  excluded 
from  permanent  sites  by  predation. 

Patterns  of  Distribution  Within  Wetlands 

Few  studies  have  examined  tadpole  distribution  patterns 
within  breeding  sites.  Diaz-Paniagua  (1987)  examined  the 
spatial  distribution  of  six  species  occurring  in  16 
temporary  ponds.  Dip-net  samples  were  taken  in  5 depth 
rones.  Vegetation  structure  varied  along  the  depth  gradient 
from  shallow  zones  characterized  by  dense  grass  cover  to 
deeper  areas  with  little  or  no  vegetation.  Along  the  depth 
gradient,  vegetation  changed  from  emergent  in  shallow  water 
to  submergent  in  deep  water.  Few  tadpoles  of  any  species 


i- vegeta ted  zones.  Although  most 


species  were  found  in  several  of  the  vegetated  zones, 
interspecific  differences  in  spatial  distribution  were 
detected.  Three  species  were  most  abundant  in  the  shallow 
zones,  two  species  were  most  abundant  in  the  emergent  zone 
and  the  remaining  species  was  most  common  in  the  submergent 

Larval  assemblages  also  may  separate  along  thermal 
(Dupre  and  Petranka  1985)  and  water  column  (Heyer  1973, 

1976)  strata.  Ontogenetic  changes  in  microhabitat 
association  have  been  found  relative  to  size  class  in  Rana 
utricularia  (Alford  and  Crump  1983) . Alford  (1986)  found 
evidence  for  both  intraspecific  (i.e.,  size  class)  and 
interspecific  spatial  partitioning  among  four  anuran  species 
within  a single  temporary  pond  site.  Two  species,  Pseudacris 
ocularis  and  Pseudacris  ornata.  separated  relative  to  cover 
and  water  temperature,  and  their  responses  to  these 
environmental  features  changed  over  time.  Alford  (1986) 
suggested  that  these  responses  were  probably  largely  due  to 
competition  (i.e.,  avoidance  relative  to  species  or  size 
class) . However,  he  did  not  examine  the  potential  influence 
ijt  predation. 

Differences  in  predator  abundance  or  composition  also 
may  exist  within  a breeding  site.  Although  poorly  studied  in 
the  tadpole-predator  system,  the  potential  for  spatial 
segregation  between  predators  and  prey  to  be  an  important 


£actor  within  wetlands  sites  was  demonstrated  by  results 
from  Banks  and  Beebee  (1988)  who  found  low  spatial  overlap 
between  Bufo  calami ta  tadpoles  and  the  invertebrate 
predators  Notonecta  sp.  and  Dvtiscus  sp.  Tadpoles  were 
concentrated  in  shallow  water,  whereas  predators  were 
collected  most  frequently  in  deeper  areas.  Because  few 
studies  have  examined  spatial  distribution  of  tadpoles 
within  breeding  sites,  it  is  not  surprising  that  there  is 
limited  information  about  the  spatial  relationships  of 
tadpole  and  predator  distributions. 

Spatial  overlap  at  the  within-wetland  scale  among 
tadpoles  and  aquatic  insect  predators  may  be  high  if  both 
groups  select  habitats  within  the  wetland  based  on  similar 
abiotic  features  such  as  water  depth,  temperature,  or  plant 
cover.  In  contrast,  differential  selection  criteria  could 
lead  to  segregation  of  tadpoles  and  predators.  For  example, 
females  of  some  odonates  require  emergent  vegetation  for  egg 
laying  (Corbet  19801 . In  this  circumstance,  tadpole  survival 
may  be  largely  influenced  by  patterns  of  microhabitat  use 
and  factors  such  as  habitat  structure.  A better 
understanding  of  how  patterns  of  tadpoles  and  predators  vary 

abiotic  and  biotic  (e.g.,  vegetation  structure)  factors, 
should  provide  valuable  insight  into  the  potential  for  such 
patterns  to  generate  differences  in  species  richness  and 
assemblage  structure  among  various  sites. 
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CHAPTER  2 

USE  OP  TEMPORARY  WETLANDS  BY  AMIRANS  IN  A MODIFIED 
LANDSCAPE:  EFFECTS  OF  HYDROLOGY  AND  WETLAND  SIZE 

Introduct ion 

Patterns  o£  distribution  of  larval  anurans  among 
aquatic  sites  are  largely  determined  by  two  related  factors: 
the  permanence  of  the  site  (i.e.,  hydrology)  and  the 
composition  of  predators  (Skelly  1995;  Werner  and  McPeek 
1990).  These  factors  are  strongly  interrelated  because 
predator  composition  changes  along  hydrologic  gradients. 

Pish  predators  dominate  in  permanent  siLes;  however,  sysLems 
that  have  an  annual  cycle  of  tilling  and  drying  prevent 
predatory  fish  populations  from  becoming  established. 

Interspecific  differences  in  life  history  trails  among 
anurans,  such  os  length  of  the  larval  developmental  period 
and  ability  to  breed  successfully  with  fish  predators,  limit 
the  range  of  wetlands  within  which  a species  can  breed 
successfully.  Species  with  long  larval  periods  are  excluded 
passively  from  siLes  with  short  hydroperiods,  whereas  most 
species  that  breed  in  temporary  sites  do  not  possess  the 
necessary  chemical  or  behavioral  characteristics  tor 
avoiding  the  heavy  predaLlon  pressure  found  in  permanent 
sites  (Kats  et  al.  1988;  Woodward  1983).  Because  of 
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whether  wetland  hydrology  or  area  aftecLed  use  by  anurans. 
further,  1 compared  breeding  activily  during  ( summer 
periods  that  varied  in  summer  rainfall  patterns. 
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Methods 

Study  Site 

Research  was  conducted  at  the 
Research  Center  (MAERC) 

Florida.  The  site  is  an 
4,800  ha  of  improved  ani 
wilh  emergent  wetlands  , 

(Sabal  palmetto!  hammock. 

The  ranch  was  ditched  extensively  in  the  1900's  to 
enhance  drainage,  and  most  wetlands  are  connected  Lo  Lhe 
vast  system  ol  ditches  that  ultimately  enters  into  Harney 
Pond  Canal,  which  flows  into  Lake  Okeechobee.  (Figure  2-11 
Most  wetlands  are  properly  characterised  as  emergent 
freshwater  marsh.  Thick-stemmed  emergents,  such  as 
Poritederia  cordata  and  Saqittarla  lanci  to  I 
species  in  deeper 
punctatum  and  Alteranthera  phlioxeroides  as  sub-dominant 
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ay  lower-stature  emergents  such  as  llacopa 


0 dominant 


•in  I 1,  i i Mil'll  i d I'd  I ! Id  did 


biodig  vil  tj  i m d i id  - 


For  the  purposes  of  this  study,  the  nine  samples  provided  a 
simple  measure  u t breeding  use  (presence/absence),  leluLive 
abundance,  and  species  richness  at  each  site. 

Wetland  area  was  estimated  from  data  previously  entered 
into  a Geographic  Information  System.  Wetland  area  varied 
f rom  0.11  to  2.21  ha  (Figure  2-1')  . Water  depth  was  taken 
from  a permanent  PVC  stake  at  the  center  of  each  wetland  at 
. 

all-weather  gage  at  the  research  center.  Ilydroperiod  is 
usually  measured  as  the  total  number  of  days  a wetland  held 
standing  water.  I took  measurements  every  three  weeks; 
therefore,  1 do  not  have  actual  dates  of  wetland  filling  and 
drying.  During  dry-down  periods,  I <ii>-i  tied  we  I lands  between 
sampling  periods  to  determine  whether  wetlands  actually 
dried  completely.  Therefore,  1 defined  ilydroperiod  as  the 
number  of  weeks  a wetland  held  standing  water. 

Statistical  Analyses 

Data  for  each  site  were  treated  as  independent  samples, 
as  were  data  from  each  summer  breeding  period.  I used 
Pearson  Product  Moment  correlation  to  determine  whether 

wetland  hydroperiod  or  size. 

During  199b,  some  wetlands  were  allected  by  spill-over 
of  water  from  diLches  containing  fish  predators  (e.g.f 
Lepomis  sp.l.  To  examine  the  effect  of  this  disturbance  on 


assemblage  composition,  I used  a chi-square  goodness  of  fit 
test  to  examine  yearly  differences  in  the  number  of  ladpoles 
that  were  from  species  that  will  or  will  not  breed  in  sites 
with  fisli  predators.  Because  not  ali  wetlands  were  altecfed 
by  this  disturbance,  I used  the  proportion  of  tadpoles  in 
each  of  the  two  categories  (i.e.,  temporary  breeding  species 
and  permanent  breeding  species)  in  unimpacted  wetlands  to 
generate  enpected  values  tor  tadpole  numbers  in  wetlands 
that  were  impacted. 


Species  Richness  and  Abundance 

A loLal  of  3,67(1  tadpoles  from  II  species  was 
collected.  Of  the  11  species  caught,  only  the  southern 
leopard  frog  (liana  utricular ia)  was  found  al  every  wetland 
(Table  2-1) . The  squirrel  treefrog  I Hyla  souirella)  was  the 
ne:;l  most  widely  distributed.  In  contrast,  the  oak  Load 
IBufo  guercicus)  and  southern  chorus  frog  (Pseudacris 
niqrita)  were  collected  from  only  two  wetlands  each. 

R.  utricularia  was  the  most  abundant  species, 
accounting  for  about  44.  of  all  captures  (Table  2-1;  figure 

grass  frog  (Pseudacris  ocularis) , pig  frog  (Kuna  qryi iol  and 
pinewoods  treefrog  (Hyla  femoral  is)  were  moderately 
abundant.  The  remaining  species  eacti  accounted  loi  Vf  oi 
less  of  the  total  sample. 


d wilh  wetland  size  (Figure  2-4;  r=0.bi 
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banks  and  spilling  water  into  adjacent  wetlands.  Comparisons 
of  anuran  breeding  activity  between  1991  and  199b  indicated 
that  species  composition  was  affected  by  this  flooding  event 
(figure  2-9) . H.  sgnireila  and  ■.  irolinen:  i . iwo  species 
that  breed  exclusively  in  temporary  sites,  were  less 
abundant  in  199b  compared  to  199-1.  In  contrast,  species  that 
associate  more  strongly  with  permanent  sites,  such  as  11. 
qrylio  and  A.  gryilus  were  more  abundant  in  199b.  Wetlands 
that  were  unimpacted  by  flooding  had  similar  proportions  of 
tadpoles  irom  species  that  will  or  will  not  breed  in  sites 
with  fish  predators  (Figure  2-101.  In  contrast,  wetlands 
that  were  impacted  by  spill-over  from  ditches  diltered 
greatly  in  composition  between  1994  and  199b  (Figure  2-10; 

=1008,  df=3,  pCO.OOOl). 


Discussion 

Changes  due  to  ditching  at  MAERC  have  resulted  in 
changes  in  wetland  hydrology,  and  in  a probable  shift  in 
assemblage  structure  and  species  composition  that  favors 
anurans  that  associate  more  strongly  wiLh  wetlands  lacking 
fish  predators.  Although  historical  data  are  not  available 

patterns  that  are  probably  similar  Lo  pre-ditching 
conditions  al  MAERC  such  as  the  Northern  Everglades,  R. 
gryl io,  A.  gryi 1 ns,  and  H.  cinerea  dominate  the  tadpole 
community  (Babbitt  et  al.,  unpubl.  data).  These  species  are 
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often  associated  with  longer  hydroperiod  sites  that  contain 
iish  predators,  and  were  probably  much  more  abundant  than 
they  are  now  (Carr  1940;  Duellman  and  Schwartz  1998) . 
Although  each  of  these  species  bred  in  some  of  the  weLlands 
I sampled,  they  concentrated  most  of  their  breeding  activity 
in  Lhe  larger  ditches  and  an  artificial  late  (pers.  obs.). 

Most  of  the  wetlands  on  site,  including  those  sampled 
during  this  study,  are  intermediate  in  hydroperiod  and  are 
best  described  as  annually  or  semi-annually  drying  temporary 
weLlands.  By  viewing  habitat  drying  as  a disturbance,  we 
would  predict  that  the  temporary  wetlands  at  MAEKC  would 
have  the  highest  species  richness  (compared  to  permanent  or 
very  ephemeral  sites)  based  on  predictions  from  the 
intermediate  disturbance  hypothesis  (sensu  Connell  1918) . 
Both  lleyer  et  al.  (1975)  and  Wilbur  (1984)  argued  that  the 
highest  species  diversity  would  be  at  intermediate  values 
along  a hydrologic  gradient.  Only  one  species  that  occurs  at 
MAERC,  the  gopher  frog  (Rana  capito) , was  not  collected 
during  the  study.  This  species  has  an  extremely  limited 
distribution  on-site,  but  it  does  breed  in  temporary 
wetlands  similar  to  those  used  in  the  sLudy  (pers.  obs.). 


wetlands  that  were  intermediate  in  hydroperiods.  In 
comparison,  only  a sub-set  oi  anuran  species  at  MAERC  breed 
in  permanent  wetlands  or  wetlands  witii  very  short 


hydroperiods. 


habitat  for  R.  utricularia.  R.  utricularia  was  the  most 
abundant  species  and  had  the  broadest  distribution.  This 
species  also  has  been  found  to  be  numerically  dominant  or 
abundant  in  censuses  at  other  temporary  weLiand  sites  in 
Florida  (Dalyrumple  1988;  Enge  and  Marion  1986;  O'Neill 
1995;  Vickers  et  al.  198b). 

Viewed  at  the  landscape  level  (i.e.,  MAERC)  the 
ditching  of  MAERC  probably  increased  anuran  species 
richness.  The  effects  of  wetland  alterations  on  anurans  al 
MAERC  differ  from  those  that  occur  at  many  other  sites 
because  this  site  was  altered  from  a large  {usually) 
permanent  aqua L i c mosaic  similar  to  the  Everglades  to  a 
terrestrial  system  containing  numerous  wetlands  of  varying 
hydroperiods.  Reductions  in  hydroperiod  aL  sites  that  are 
already  temporary  can  result  in  decreased  species  richness 
or  abundance  because  Lruncated  hydroperiods  preclude  some 
species  from  successfully  metamorphosing  (Rechmann  et  al. 
1989;  Rowe  and  bunson  1995;  Semiilsch  1987;  Wilbur  1987). 

The  effects  of  ditching  of  wetlands  on  anuran  species 
richness  and  breeding  success  depend  on  both  the  original 
and  resultant  conditions  of  a site. 

abundance  was  related  to  wetland  hydroperiod.  These  results 
contrast  with  those  of  Rechmann  el  al.  11989)  who  found  that 
hydroperiod  was  positively  correlated  with  species  richness 
and  the  abundance  of  metamorphosing  amphibians,  in  their 


study,  hydroperiods  varied  almost  by  a factor  of  five  (i.e., 
•,!{  Lo  263  days),  liydroperiods  in  the  current  study  varied  by 
only  a factor  of  1.4,  from  49  to  69  weeks,  further,  the 
hydroperiod  of  these  wetlands  was  much  longer  Ilian  the 
minimum  required  for  metamorphosis  of  all  Lhe  anuran  species 
occurring  at  (MERC  (Ugas  19611;  Will..,.  191.7) . In  a live-year 
study,  Dodd  (1992)  did  not  find  a relationship  between 
habitat  duration  and  species  richness  at  a single  temporary 
wetland  in  north-central  Florida.  That  study  was  conducted 
during  a drought,  and  visitation  by  aduiL  anurans  was  not 
necessarily  associated  with  breeding  activity.  Large 
differences  in  hydrology  among  sites  (or  among  yearsl  may  be 
required  to  generate  significant  differences  in  species 
richness.  Relevant  differences  may  include  sites  with 
hydroperiods  that  are  too  short  for  some  species  to 
successfully  produce  metamorphs,  and  permanent  sites  with 
fish  predators.  It  is  likely  LI, at  the  wetlands  in  the 
current  study  provided  similar  lwibil.,1  sui Labi  1 ,ly,  .elal.ve 
to  hydrology,  and  that  any  variation  in  use  patterns  by 
adults  was  largely  related  to  other  factors. 

Species  richness  was  positively  related  Lo  wetland 
sir.e.  Results  from  studies  till. I have  onnm.ned  I In* 
relationship  between  wetland  size  and  species  richness  or 
abundance  are  equivocal.  Whereas  some  researchers  have  lound 
.1  positive  relationship  between  wetland  area  and  species 
richness  (Kutka  and  Bachmann  19911;  Laan  and  Verboom  1990), 
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others  have  found  no  relation  aL  all  (Dias-l’aniagua  1990; 
Laan  and  Verboom  1990;  Richter  and  Azous  19991 . Again, 
differences  in  these  studies  may  represent  differences  in 
scale,  as  well  as  differences  in  oLher  characteristics  tliat 
influence  use  by  adults  (l.ann  and  Verboom  1990;  Richter  and 
Azous  1999).  because  species  richness  aL  a particular  site 
is  a community  metric  and  therefore  a product  of  the  varying 
responses  of  community  members,  iL  may  not  be  surprising 
Lhat  the  above  studies  provide  differing  results. 

Within  a particular  class  ol  weLiands,  such  as  the 
temporary  wetlands  with  relatively  long  hydroperiods  sampled 
in  this  study,  year-to-year  changes  in  hydrology  driven  by 
meteorological  conditions  may  have  a larger  influence  on 
breeding  use  than  many  other  factors.  The  largest  influence 
on  anuran  use  of  wetland  sites  at  MAKKC  was  amount  and 
pattern  of  rainfall.  During  Lhe  summer  drought  conditions  ol 
1993,  no  breeding  took  place  in  any  temporary  wetlands 
because  they  were  dry.  In  1994,  essentially  all  non- 
permanent sites  provided  suitable  habitat  for  species  that 
breed  in  temporary  weLiands  because  all  sites  dried  down 
prior  to  the  summer  rains.  This  eliminated  any  fish 
predators  that  would  have  made  these  sites  unsuitable  for  a 
majority  of  species.  These  two  years  contrasted  further  with 
1995,  when  some  wetlands  dried  completely  before  lhe  summer- 
rains  but  others  did  not  completely  dry.  Further,  over-spill 
from  adjacent  fish-containing  ditches  apparently  caused  some 


anurans  to  avoid  using  sites  thal  otherwise  would  have  been 
suitable. 

Although  the  mechanism  through  which  adult  amphibians 
are  able  to  identify  sites  that  contain  fish  is  not  known, 
experimental  choice  studies  have  demonstrated  that  adults 
are  capable  of  differentiating  between  sites  that  do  and  do 
not  contain  fish  predators  IHopey  and  Petranka  1994; 
KesitariLs  and  Wilbur  1991).  differences  in  species 
composition  among  siLes  with  and  without  fish  predators  are 
due  (at  least  partly)  Lo  adulL  selection,  not  simply 
predation  on  the  tadpoles  of  certain  species. 

The  impacts  of  over-spill  of  waLer  from  fish-containing 
ditches  into  otherwise  isolated  wetlands  are  transient 
because  wetlands  usually  dry  on  an  annual  basis.  Further, 
alternative  sites  are  available  for  anurans  that  avoid 
breeding  in  siLes  with  predatory  lish.  For  example,  although 
species  such  as  H.  squi rel la  and  C.  carol i nons is  avoided 
breeding  in  wetlands  affected  by  over-spill  Iron  fish- 
containing  ditches,  these  species  bred  in  adjacenL  sites 
that  were  unaffected  (pers.  obs.). 

In  contrast,  low  rainfall  can  result  in  missed  breeding 

exacerbates  the  dry  conditions.  During  seasonal  droughts, 
the  only  species  that  can  successiuily  breed  are  those  that 
breed  at  permanent  sites  with  fish  predators.  Whereas 
periodic  seasonal  droughts  may  have  only  a small  iinpacL  on 


population  numbers,  the  effects  of  long-term  drought  may  be 
more  profound.  Few  studies  have  examined  amphibian 
populations  during  drought,  and  the  comple::  and  erratic 
patterns  of  population  number  of  amphibians  make  assessment 
of  the  impacts  of  drought  difficult  (Dodd  1992) . However, 
some  species  appear  to  be  capable  of  surviving  long-term 
droughts  (Dodd  1995) . Evidence  that  at  least  some 
individuals  within  a population  may  be  opportunistic  in 
selecting  breeding  sites,  rather  than  philopatric,  suggests 
that  this  may  be  an  important  mechanism  for  maintaining  a 
breeding  population  during  drought  conditions  (Dodd  1995) . 
Opportunistic  use  of  breeding  sites  also  appeared  to  be 
important  in  areas  affected  by  flooding  during  1995  (pers. 
obs.).  Particularly  in  landscapes  such  as  Florida  where 
wetlands  are  numerous,  opportunistic  use  of  breeding  sites 
may  be  more  important  than  previously  recognized  (also  see 
Dodd  1995) . 

These  data  suggest  that  temporary  wetlands  on  this  site 
provide  dynamic  habitats  that  offer  varying  breeding 
opportunities  and  larval  developmental  conditions  that  are 
highly  dependent  on  meteorological  conditions.  The  effect 
ol  lliis  is  a spa Lialiy- temporally  dynamic  system  resulting 
in  differing  assemblage  structure  and  composition  at 
particular  sites.  Because  MAERC  contains  numerous  wetlands 
and  ditches  that  vary  in  hydrology  from  ephemeral  to 
permanent,  most  species  can  find  suitable  breeding  habitat 
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except  under  extreme  drought  conditions.  During  such 
conditions,  ditching  exacerbates  the  drought  conditions  by 
increasing  drainage,  lowering  the  water  table,  and 
prolonging  dry  conditions. 

In  conclusion,  historical  changes  to  the  landscape  have 
probably  increased  local  anuran  species  richness  and  altered 
tadpole  assemblage  structure  at  the  wetlands  at  MAERC.  An 
important  aspect  that  was  beyond  the  scope  o£  this  study  is 
the  effects  of  past  and  current  land  management  on  the 
amount  and  distribution  of  varying  terrestrial  habitats.  The 
surrounding  upland  matrix  can  have  a large  influence  on 
wetland  use  as  breeding  sites,  particularly  for  amphibians 
that  have  a terrestrial  adult  phase  (Kutka  and  Bachmann 
1990;  I,aan  and  Verboom  1990)  . For  example,  arboreal  species 
such  as  H.  femoralis  appeared  to  be  locally  abundant,  but 
limited  in  distribution.  Temporary  wetlands  near  forested 
hammocks  consistently  contained  H.  femoralis  tadpoles; 
however,  wetlands  away  from  hammocks  in  large  open  pastures 
did  not.  Consideration  of  upland  characteristics  is  an 
important  aspect  of  overall  conservation  of  anurans. 

Research  on  the  relationship  between  upland  habitat  patch 
characteristics  and  wetland  use  would  increase  our 
understanding  of  breeding  use  patterns,  as  well  as  our 
ability  to  provide  guidance  for  managing  anuran  populations. 
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Figure  2-1.  Map  of  the  Macflrthur  Agro-Ecology  Research 
Center  showing  the  system  of  ditches  that  enhance  drainage. 
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Figure  2-4.  Relationship  between  wetland  size  and  anuran 
species  richness  (top)  and  abundance  (boLlom) . Abundance  and 
richness  based  on  total  captures  from  all  0.25m'  traps 
within  each  weLland. 
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Figure  .1-5.  Relationship  between  wetland  hydroperiod  anc 
anuran  species  richness  (top)  and  abundance  (bottom). 
Abundance  and  richness  based  on  total  captures  from  all 
0.25m'  traps  within  each  wetland. 
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Figure  2-10.  Comparison  of 

the  1994  and  1995  sv 
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sample  composition  betwee 
breeding  seasons  at  wetlands  that 
were  m=oi  <mu  .•»-  .n=5)  impacted  by  flooding  from 

ditches  containing  fish  predators.  Species  categorized  as 
temporary  site  breeders  were  species  that  do  not  breed  in 
sites  with  fish  predators.  Rana  grylio,  Rana  utricularia, 
Acris  qrvllus,  and  Hyla  cinerea  will  breed  with  fish 
predators  and  were  categorized  as  permanent  s 


e breeders. 
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CHAPTER  3 

AMD  TEMPORAL  DYNAMICS  OF  TADPOLES  AND  AQUATIC  INSECT 
PREDATORS  DEVELOPING  IN  TEMPORARY  WETLANDS 


Introduction 

Predation  is  a major  source  of  mortality  in  anuran 
tadpoles  (Calef  1973;  Heyer  et  al.  1975;  Smith  1983),  and  is 
thought  to  have  a large  influence  on  breeding  site  selection 
by  adult  anurans  (Magnusson  and  Hero  1991) . Fish  are  major 
predators  of  tadpoles  in  permanent  aquatic  sites,  and  many 
species  appear  to  avoid  breeding  in  sites  that  contain  fish 
predators  (Hopey  and  Petranka  1994;  Resetarits  and  Wilbur 
1991) . Although  ephemeral  wetland  sites  generally  do  not 
contain  fish  predators,  the  short  duration  of  these  sites 
makes  them  unsuitable  for  species  with  long  larval 
developmental  periods  IWilbur  1980,  1984,  1987) . Thus 
differences  in  habitat  duration  and  predator  composition, 
coupled  with  interspecific  differences  in  the  selection  of 
breeding  sites  by  adults,  can  generate  differences  in 
species  composition  of  larval  anurans  along  gradients  of 
habitat  permanence  (Werner  and  McPeek  1994). 

However,  factors  that  influence  tadpole  success,  such 
as  abundance  and  composition  of  competitors  or  predators, 
vary  in  space  not  only  among  wetland  sites  (Gascon  1992; 
Smith  1983;  Werner  and  McPeek  1994;  Woodward  1983)  but  also 
40 


among  patches  within  wetlands  (Banks  and  Beebee  1988;  Diehl 
1988,  1992) . Therefore,  the  distribution  patterns  of 
tadpoles  and  predators  at  both  the  macrohabitat  (breeding 
site)  and  microhabitat  (patch  site)  scale  can  have  a large 
influence  on  tadpole  growth  and  survival,  and  thereby 
assemblage  structure. 

Although  the  importance  of  considering  scale  has  been 
widely  discussed  (Wiens  1989,  1992),  few  studies  of  tadpole 
ecology  have  addressed  this  issue.  Processes  that  are 
important  in  population  regulation  at  one  scale  of 
investigation  may  fail  to  operate  at  other  scales  (Wiens 
1992) . For  example,  whereas  predator  composition  has  a large 
influence  on  anuran  breeding  site  selection,  the  spatial 
dynamics  of  predators  and  tadpoles  at  the  within-wetland 
scale  are  poorly  understood.  If  tadpoles  occupy  patches 
that  provide  increased  protection  from  predators,  such  as 
areas  with  high  cover  (Babbitt  and  Jordan  1996;  Chapters  4 
and  6),  or  areas  with  lower  predator  densities  (Banks  and 
Beebee  1988),  then  differential  habitat  use  within  sites  may 
be  an  important  feature  influencing  predation  rates,  and 
possibly  assemblage  structure  (Morin  1986) . 

Spatial  overlap  at  the  within-wetland  scale  among 
tadpoles  and  predators  may  be  high  if  both  groups  select 
microhabitats  based  on  similar  abiotic  features  such  as 
water  depth,  temperature,  or  plant  cover.  Understanding 
such  patterns,  and  how  they  vary  in  space  and  time,  should 


provide  valuable  insight  into  the  potential  for  differences 
in  microhabitat  selection  to  generate  differences  in  species 
richness  and  success  among  various  sites.  For  example, 
occupation  of  vegetated  shallow  areas  of  permanent  aquatic 
sites  that  have  deeper,  poorly  vegetated  open-water  areas 
may  provide  increased  protection  from  fish  predators  (Diehl 
1992;  Mclvor  and  Odum  1988;  Werner  et  al.  1983).  How 
distribution  patterns  influence  predation  rates  on  tadpoles 
developing  in  temporary  wetlands  where  aquatic  insects  are 
Ihe  major  predators  is  largely  unknown. 

Several  studies  have  found  interspecific  differences  in 
tadpole  distribution  within  wetland  sites  (e.g,  Alford  1986; 
Diaz-Paniagua  1982;  Heyer  1976);  however,  few  studies  have 
related  tadpole  distribution  to  predator  distribution  (Banks 
and  Beebee  1988) . The  lack  of  research  on  within-site 
distribution  patterns  of  tadpoles  and  aquatic  insect 
predators  leaves  open  the  question  of  whether  differential 
distribution  patterns  within  wetlands  is  a potential 
mechanism  for  decreasing  predation  pressure  on  tadpoles. 

The  goal  of  this  study  was  to  determine  the  spatial  and 
temporal  patterns  of  anuran  larval  and  aquatic  insect 
predators  within  several  temporary  wetland  sites.  I asked 
whether  tadpoles  or  predators  were  differentially 
distributed  in  three  rones  based  on  water  depth.  Further,  I 
asked  whether  habitat  features  influenced  tadpole 
distribution  to  determine  whether  tadpoles  differentially 
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occupied  patches  based  on  physical  or  chemical  habitat 
features.  By  examining  distribution  patterns  within 
wetlands,  I was  able  to  determine  whether  differences  in 
distribution  could  function  as  a mechanism,  in  addition  t 
adult  breeding  site  selection,  for  decreasing  predation 
pressure  on  larval  anurans. 


Study  Site 

Research  was  conducted  at  the  MacArthur  Agro-Ecology 
Research  Center  in  Lake  Placid,  Highlands  County,  Florida. 
The  site  is  an  active  cattle  ranch  consisting  of  4,800  ha  o 
improved  and  semi-improved  pasture  interspersed  with 
emergent  wetlands  and  oak  IQuercus  virqiniana)  and  palm 
I Saba 1 palmetto)  hammock. 

The  ranch  was  ditched  extensively  in  the  196Q'S  to 
enhance  drainage,  and  most  wetlands  are  connected  to  the 
vast  system  of  ditches  that  ultimately  enters  i 
Pond  Canal,  which  flows  into  Lake  Okeechobee.  M 
are  properly  characterized  as  emergent  freshwater  marsh. 
Thick-stemmed  emergents  such  as  Pontederia  cordata  and 
Saoittaria  lancitoiia,  are  dominant  species  in  deeper  water 
with  Panicum  hemitomon.  Polygonum  punctatum  and  Alteranthera 
philoxeroides  as  sub-dominant  species.  Outer,  shallower 
areas  of  temporary  wetlands  are  characterized  by  lower- 
stature  plants  such  as  Bacopa  caroliniana,  Hydrochloa 


t wetlands 


carolinensis  and  Diodia  virqiniana. 


Sampling  Methods 

From  May  1994  through  September  1995,  I sampled  12 
wetlands  every  three  weeks  (n=23  sampling  periods) . One 
wetland  was  dropped  during  the  last  six  sampling  periods  due 
to  damage  from  cattle  trampling.  At  each  site,  I collected 
9 samples  using  a stratified  random  sampling  protocol.  Three 
replicate  samples  were  collected  in  each  of  3 concentric 
belts  corresponding  approximately  to  one  third  the  radius  of 
the  wetland. 

Tadpoles  and  macroinvertebrales  were  collected  with  a 
1/4-m'1  open-ended,  box  trap.  Box  traps  have  been  used  to 
sample  small  fish  and  macroinvertebrates  (Chick  et  al.  1992; 
Freeman  et  al.  1984;  Kushlan  1981)  and  anuran  larvae 
(Caldwell  et  al.  1980;  Calef  1973;  Pfenning  1990),  and  are 
effective  for  sampling  in  areas  with  complex  vegetation 
structure.  A bar  seine  fitting  the  diameter  of  the  box  trap 
and  fitted  with  fiberglass  insect  screening  was  swept 
through  the  water  column  (after  vegetation  had  been  cleared) 
until  3 consecutive  sweeps  yielded  no  additional  captures. 
Samples  were  preserved  in  the  field  with  10%  buffered 
formalin.  Macroinvertebrates  were  removed  from  formalin, 
rinsed  thoroughly  in  water  and  stored  in  alcohol  in  the  lab. 
Tadpoles  and  macroinvertebrates  were  identified  to  species. 
For  the  purposes  of  this  study,  the  macroinverlebrate 
assemblage  was  analyzed  as  a whole,  rather  than  by  species. 
The  assemblage  consisted  of  11  species  of  odonates  from  the 


families  fteshnidae,  Libellulidae  and  three  species  of 
Hemiptera  from  the  families  Belostomatidae  and  Naucoridae. 

Water  depth  was  measured  to  the  nearest  centimeter 
inside  each  trap  prior  to  removal  of  plants.  Foliar  coverage 
of  above-ground  vegetation  was  estimated  ocularly,  and  then 
plants  were  removed  and  shaken  to  dislodge  any  macrofauna. 
Grass  clippers  were  used  to  clip  plants  as  necessary  to 
remove  all  vegetation.  Vegetation  was  placed  in  a mesh  bag, 
which  was  spun  in  the  air  20  times  to  remove  excess  water 
prior  to  measuring  biomass  to  the  nearest  0.1  kg.  Water 
temperature  ("G>  and  pH  measurements  were  taken  adjacent  to 
each  trap  at  0.1  m below  the  water  surface. 

Statistical  Analyses 

Split-plot  analysis  of  variance  (ANOVA)  was  used  to 
test  for  effects  of  wetland  site,  sampling  period  (time), 
wetland  rone  and  the  wetland  ::  zone  and  sampling  period  :: 
zone  interactions  on  the  habitat  features  water  depth,  water 
temperature,  pH,  cover,  and  plant  biomass,  and  on  tadpole 
and  macroinvertebrate  density  and  biomass.  In  addition, 
separate  analyses  were  done  on  the  three  most  abundant 
anuran  species.  Values  from  the  three  samples  taken  within 
each  zone  were  averaged  to  calculate  a mean  value  for  each 
zone  for  all  parameters.  Tadpole  and  macroinvertebrate 
density  and  biomass  values  were  log  transformed  (log  (k+11) 
prior  to  analysis  to  reduce  skewness.  Examination  of 
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residuals  for  the  measures  of  environmental  parameters 
indicated  that  transformations  were  not  necessary  (Sokal  and 
Rolf  1981).  The  main  effects  of  wetland  and  zone  were  tested 
over  the  split-plot  error  (i.e.,  wetland  x zone),  whereas 
the  time  and  time  x zone  effects  were  tested  over  the  whole 
plot  error  (i.e.,  mean  square  error).  The  Bonferonni-Dunn 
post  hoc  multiple  comparisons  tests  was  used  to  test  for 
differences  among  the  three  zones  when  a significant 
(p<0 . 05)  zone  effect  was  detected  (Day  and  Quinn  1989). 

Correlation  analyses  were  performed  to  examine  the 
relationships  between  the  abundance  of  all  tadpoles 
combined,  as  well  as  the  abundance  of  the  three  most  common 
anurans  (H.  squirella.  R.  utricularia,  and  8.  ocularis) , 
with  habitat  variables  and  macroinvertebrate  abundance  and 
biomass  (Sokal  and  Rolf  1981) . 


Results 

Physical  Parameters 

Water  depth,  water  LemperaLure,  pH,  plant  cover,  and 
plant  biomass  all  varied  significantly  among  wetlands  and 
across  time  (Table  3-1;  Figures  3-1  through  3-5) . The  main 
etfecL  of  Lime,  which  reflecLs  seasonal  differences, 
accounted  for  most  of  the  variation  in  water  depth  (37%)  and 
water  temperature  (58%).  In  contrast,  wetland  explained  a 
majority  of  the  variation  in  pH  (73%),  plant  cover  (36%), 
and  plant  biomass  (29%).  Wetlands  that  were  adjacent  to 
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upland  hammocks  tended  to  have  lower  pH  values  compared  to 
wetlands  that  were  surrounded  by  pasture.  Wetland 
differences  in  plant  cover  and  biomass  indicate  differences 
in  plant  communities,  as  well  as  differences  in  the  areal 
extent  of  coverage  of  thicker-stemmed  emergents. 

Significant  differences  among  zones  were  detected  for 
water  depth,  pH,  plant  cover  and  plant  biomass,  but  not  for 
water  temperature  (Figures  3-1  through  3-5).  Both  water 
depth  and  pH  increased  progressively  from  the  outer  to  the 
inner  zone.  Whereas  the  zone  effect  explained  22%  of  the 
variation  in  water  depth,  it  accounted  for  only  1%  of  the 
variation  in  pH. 

Cover  levels  were  higher  in  the  outer  wetland  zones 
compared  to  the  middle  and  inner  zones  (Figure  3-4) . Cover 
in  the  outer  zone  was  nearly  100%  at  most  wetlands 
throughout  the  study.  Plant  biomass  also  differed  among 
zones;  however,  in  contrast  to  cover  where  highest  levels 
were  in  the  outer  zone,  plant  biomass  was  highest  in  the 
inner  zone  (Figure  3-5) . Whereas  difference  among  zones 
accounted  for  20%  of  the  variation  in  cover, 
for  only  5%  of  the  variation  in  biomass. 

Few  interactions  were  significant,  and  m 
explained  very  little  of  the  variation  in  the 
The  significant  zone  x wetland  interaction  for  plant  a 
and  biomass  was  due  to  differences  among  wetlands  in  t 
extent  of  the  development  of  the  inner  thick-stemmed 


L accounted 


L that  were 


emergent 


Tadpole  and  Predator  Distribution 

Tadpole  densities  varied  significantly  among  wetland, 
time,  and  wetland  zone  (Table  3-2,  Figure  3-6) . Wetland  and 
time  explained  9%  and  281  of  the  variation  in  tadpole 
density,  respectively.  Densities  were  higher  in  the  outer 
zone  compared  to  both  the  inner  and  middle  zones,  but  did 
not  differ  significantly  between  the  inner  and  middle  zones 
(Figure  3-6) . However,  zone  explained  less  than  one  percent 
of  the  variation  in  tadpole  densities.  Such  differences  are 
probably  biologically  insignificant.  In  contrast  to  density, 
tadpole  biomass  did  not  differ  significantly  among  zones 
(Table  3-2,  Figure  3-7).  Differences  among  wetlands 
explained  91  of  the  variation  in  tadpole  biomass  and  time 
explained  231.  Although  tadpole  numbers  were  highest  at  the 
beginning  of  each  summer  breeding  period,  tadpole  biomass 
was  highest  during  the  winter  months.  This  difference 
reflects  the  dominance  of  the  relatively  larger  tadpoles  of 
R.  utricularia  in  winter  samples. 

The  pattern  for  macroinvertebrate  predators  differed 
from  that  of  tadpoles.  Whereas  zone  differences  were  found 
for  tadpole  density  but  not  biomass,  the  opposite  was  true 
lor  predaLois  (Table  3-2).  MacroinverCebrale  numbers  did  nol 
differ  among  zones;  however,  macroinvertebrate  biomass  was 
higher  in  the  middle  zone  compared  to  the  outer  zone. 
(Figures  3-8  and  3-9).  Again,  although  zone  differences  in 
biomass  were  significant,  they  explained  less  than  one 
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percent  of  the  variation  in  biomass.  In  contrast  to  tadpole 
density  and  biomass  where  time  explained  most  of  the 
variation,  wetland  explained  a higher  percentage  of  the 
variation  in  macroinvertebrate  number  and  biomass  than  did 
Lime.  Wetland  accounted  for  34%  and  24%  of  the  variation  in 
macroinvertebrate  density  and  biomass,  respectively,  whereas 
time  explained  11%  and  17%.  Although  the  numbers  of 
macroinvertebrates  varied  significantly  among  wetlands, 
wetland  size  was  not  correlated  with  predator  numbers 
(R=0.258,  p-0.441. 

Species-Specific  Responses 

Three  species,  H.  squirella,  R.  utricularia,  and  P. 
ocularis,  were  abundant  enough  to  examine  separately. 

Wetland  and  time  had  sfqni Ficant  effects  on  density  of  each 
species  (Table  3-3;  Figures  3-10  through  3-12) . For  H. 
souirella,  time  accounted  for  a majority  of  variation  in 
tadpole  numbers  (40%) . Wetland  explained  only  4%  of  the 
variation  in  H.  souirella  numbers.  The  main  effects  of 
wetland  and  time  explained  18%  and  16%,  respectively,  of  the 
variation  in  R.  utricularia  densities.  These  effects  each 
explained  12%  of  Lhe  variation  in  P.  ocularis  densities.  R. 
utricularia  and  P.  ocularis  breed  throughout  the  year,  so  it 
is  not  surprising  that  time  explained  less  of  the  variation 
in  density  of  these  two  species  compared  to  H.  squirella, 
which  is  strictly  a summer 


breeder.  Differences  in  the 
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pattern  of  distribution  among  zones  were  not  significant  for 
any  of  these  species. 

Relationship  to  Physical  Parameters 

Although  significant  relationships  were  found  between 
tadpole  numbers  and  many  environmental  variables,  these 
relationships  were  very  weak  (Table  3-4) . The  strongest 
relationship  was  a negaLive  correlation  with  waLer  depth. 
More  complex  patterns  emerge,  and  some  stronger 
relationships  were  detected,  when  the  three  most  abundant 
species  are  examined  separately.  H.  squirella  had  a strong 
negative  relationship  with  water  depth  and  a strong  positive 
relationship  with  temperature  and  pH  (Table  3-4) . In 
addition,  H.  squirella  was  negatively  correlated  with  plant 
biomass.  R.  utricularia  also  had  a negative  relationship 
with  water  depth  but  this  relationship  was  weak.  In  contrast 
to  H.  squirella,  R.  utricularia  was  negatively  correlated 
with  temperature.  The  differences  between  these  two  species 
reflect  the  interspecific  differences  in  seasonality  of 
breeding  activity.  R.  utricularia  also  was  positively 
correlated  with  plant  cover.  P.  ocularis  numbers  were  not 
correlated  significantly  with  any  ol  the  environmental 
variables  measured. 


Relationship  to  Macroinvertebrates 


Tadpole  numbers  and  macroinvertebrate  numbers  were  not 
significantly  related,  and  only  a weak  negative  relationship 
was  found  between  tadpole  numbers  and  macroinvertebrate 
biomass  (Table  3-4) . No  significant  relationships  were  found 
between  macroinvertebrate  number  or  biomass  and  any  of  the 
three  species  examined  separately. 

Discussion 

Patterns  of  tadpole  abundance  within  wetlands  provided 
little  evidence  that  differential  occupation  within  wetland 
zones  was  an  important  anti-predator  mechanism.  Although 
differences  in  tadpole  abundance  among  wetland  zones  were 
statistically  significant,  these  differences  explained  only 
one  percent  of  the  total  variation  in  tadpole  numbers. 
Examination  of  the  entire  tadpole  assemblage  as  a group  may 
obscure  some  species-specific  patterns  that  may  exist; 
however,  separate  analyses  of  the  three  most  abundant 
species  did  not  reveal  any  differences  in  distribution 
relative  to  wetland  zone.  Because  predator  numbers  and 
biomass  also  were  distributed  relatively  equally  among 
zones,  differential  use  by  tadpoles  presumably  would  not 
have  been  an  effective  mechanism  for  avoiding  predators. 

In  addition  to  the  lack  of  differences  in  distribution 
based  on  wetland  zone,  tadpoles  showed  weak  relationships 
with  habitat  characteristics.  Examination  of  the  entire 


tadpole  assemblage  together  revealed  non-significant  or  weak 
correlations.  However,  patterns  among  the  three  species 
examined  separately  indicated  that  the  distribution  of  at 
least  one  species,  H.  squirella,  was  related  to  some 
environmental  parameters.  This  species  was  negatively 
correlated  with  depth  and  plant  biomass,  and  positively 
correlated  with  water  temperature  and  pH.  The  positive 
correlation  with  temperature  was  due  to  exclusive  summer 
breeding  rather  than  differences  within  wetlands  (which  were 
not  significant) . Further,  the  posiLive  correlation  with  pH 
was  due  to  differences  among  wetland  sites,  rather  than 
differences  within  wetlands.  Both  water  depth  and  plant 
biomass  were  significantly  lower  in  outer  wetland  rones, 
suggesting  that  the  negative  correlation  between  these 
parameters  and  H.  squirella  numbers  was  due  to  use  of  outer 
zones.  However,  ANOVA  Indicated  that  H.  squirella  numbers 
were  not  related  to  rone.  Because  H.  squirella  breeds 
shortly  after  wetland  fills  from  rainfall,  use  of  shallow 
wetlands  and  wetlands  with  less-developed  rones  of  thick- 
stemmed emergents  probably  contributed  to  the  negative 
correlations.  Thus,  correlations  of  H.  squirella  with 
habitat  features  appeared  to  be  largely  related  to  eatly 
summer  breeding  and  adult  selection  of  wetland  breeding 

R.  utricularia  numbers  were  related  only  to  three 
variables,  cover  (positive),  water 


depth  (negative)  and 
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temperature  (negative) . Similar  to  H.  squirella,  differences 
may  largely  be  related  to  differences  among  wetland  sites, 
since  R.  utricularia  numbers  also  did  not  differ 
significantly  among  zones.  The  negative  correlation  with 
temperature  reflects  winter  breeding  activity.  Finally,  P. 
ocularis  numbers  were  not  related  to  any  variables  measured. 

The  results  suggest  that  within  a class  of  wetland, 
(i.e.,  temporary),  predator  avoidance  may  not  be  related  to 
microhabitat  selection,  at  least  at  the  scale  examined  in 
this  study.  At  the  scale  of  wetland  site,  Gascon  (1991) 
found  that  species-specific  responses  to  habitat  features  at 
several  temporary  wetland  sites  were  highly  variable  and 
were  not  consistent  from  year  to  year.  Thus,  whereas  species 
may  breed  only  in  temporary  sites,  they  may  use  sites  that 
are  poorly  vegetated  or  well  vegetated  (see  also  Heyer 
1976) . 

Some  studies  that  have  examined  distribution  patterns 
of  tadpoles  have  found  differential  distribution  within 
sites.  Generally,  the  patterns  observed  have  been  higher 
occupation  of  vegetated  areas  and  lower  occupation  of  deeper 
open-water  areas  or  areas  with  very  low  cover  (Alford  1986; 

Loschenkohl  1986) . Because  wetland  sites  in  the  current 
study  were  vegetated  throughout,  differences  in  results  with 
previous  studies  may  indicate  that  large  differences  in 
habitat  features,  such  as  cover  or  no  cover,  may  be 


necessary  tor  generating  differential  distributions. 
Although  cover  levels  did  vary  significantly  among  zones, 
these  differences  may  not  have  been  large  enough  to  affect 
microhabitat  use. 

Observed  patterns  of  tadpole  distribution  may  be 
generated  by  a variety  of  processes,  including  differential 
predation,  tadpole  selection  of  microhabitats,  and  adult 
selection  of  oviposition  sites.  The  relative  role  of  these, 
or  alternative  mechanisms,  in  generating  observed  patterns 
is  unclear.  Unfortunately  few  studies  that  have  examined 
tadpole  distribution  have  also  examined  predator 
distribution.  An  exception  is  work  by  Banks  and  Beebee 
(1988)  who  found  that  both  predatory  dytiscid  beetle  larvae 
and  Bufo  calami ta  tadpoles  were  non-randomly  distributed 
relative  to  water  depth  and  vegetation,  with  tadpoles 
concentrating  in  shallow  water  and  beetle  larvae  being  more 


common  in  non-vegetated,  deeper  water  areas. 

An  aspect  of  microhabitat  use  not  addressed  in  this 
study  is  differentiation  based  on  water  column  strata. 
Interspecific  differences  in  use  of  the  water  column  have 
been  observed  or  inferred  from  tadpole  morphology  (e.g.. 


1980).  Such  interspecific  differences  often  have  been 
interpreted  as  a mechanism  for  avoiding  exploitation 
competition;  however,  regardless  of  the  mechanism  generating 
this  use  pattern,  such  differences  could  result  in 


differential  predation  pressure.  For  example,  although 
aeshnid  odonates  such  as  Anax  junius  cling  to  vegetation, 
many  libellulid  odonates  are  benthic  (Pritchard  1965).  Thus, 
the  foraging  strategies  of  the  macroinvertebrates  could 
generate  interspecific  differences  in  predation  on  tadpoles 
in  patches  with  differing  predator  composition. 

Many  odonates  and  anurans  oviposit  shortly  after 
wetlands  fill  with  water;  therefore,  tadpoles  developing  in 
varying  parts  of  a wetland  may  face  similar  levels  of 
predation  pressure.  For  these  species,  the  susceptibility  to 
predators  will  largely  be  a function  of  growth  rates. 
Tadpoles  that  grow  rapidly  can  effectively  eliminate 
predation  threats  from  all  but  the  largest  predators  (Crump 
1984;  Richards  and  Bull  1990;  Travis  et  al.  1985a;  but  see 
Crump  and  Vaira  1991) . Such  a mechanism  is  probably 
important  for  species  such  as  H.  squirella.  Thus,  aithough 
occupation  of  areas  with  high  cover  may  decrease  predation 
rates  on  H.  squirella  tadpoles  (Babbitt  and  Jordan  1996; 
Chapter  6) , the  timing  of  breeding,  avoidance  of  sites  with 


for  successful  tadpole  development  in  this  species. 

Although  breeding  early  may  decrease  predation  pressure 
on  some  species,  and  may  explain  why  the  tadpoles  of  these 
species  do  not  show  strong  microhabitat  use,  R.  utricularia, 
which  is  mainly  a winter  breeder,  did  not  use  this  strategy. 
It  bred  in  sites  that  contained  both  abundant  number  of 
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predators  as  well  as  large  predators  (Figures  3-8  and  3-11; 
pets,  obs) . However,  like  H.  souirella,  this  species  did  not 
show  strong  patterns  of  microhabitat  use.  Thus,  lack  of 
microhabitat  selection  to  avoid  predation  was  not  limited  to 
early  breeding  species.  Therefore,  early  breeding  offers 
only  a partial  explanation  for  the  lack  of  microhabitat 
selection.  More  likely,  the  broad  and  even  distribution  of 
the  predators  may  make  such  a strategy  ineffective  in 
wetlands  that  do  not  vary  greatly  in  habitat  features  such 
as  plant  cover. 

Although  tadpoles  did  not  exhibit  strong  correlation 


occupation  of  patches  that  differ  does  not  affect  tadpole 
growth  and  survival  (Holomuzki  1986a;  Sredl  and  Collins 
1992;  Travis  and  Trexler  1986) . Experimental  studies  have 
demonstrated  the  importance  of  increased  cover  in  mediating 
predator-prey  interactions  in  aquatic  systems  (reviewed  in 
Heck  and  Crowder  1991),  including  interactions  beLween 
tadpoles  and  aquatic  insect  predators  (Babbitt  and  Jordan 
1996;  Chapters  4 and  6) . 

The  lack  of  strong  responses  to  habitat  features  within 
wetlands  contrasts  with  the  large  Change  in  composition 
observed  in  responses  to  the  introduction  of  water  from 
ditches  containing  fish  (Chapter  2).  The  shift  from  an 
assemblage  dominated  by  species  that  breed  in  temporary 
sites  to  one  dominated  by  species  that  will  breed  in  sites 
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containing  fish  predators  was  interpreted  as  being  due  to 
adult  breeding  site  selection  (Hopey  and  Petranka  1994; 
Resetarits  and  Wilbur  19911.  Thus,  predator  composition  is 
important  for  determining  breeding  use  at  the  macrohabitat 
level.  Although  abiotic  habitat  characteristics  may 
determine  use  by  some  species,  features  within  sites  that 
fall  within  a "suitable"  range  for  breeding  may  not  be 
important  in  determining  species  use  (Gascon  1991).  In 
addition,  use  patterns  within  sites  may  not  be  strongly 
influenced  by  predators,  except  in  sites  with  clear 
partitioning  of  habitat  features  (i.e.,  cover  versus  not 
cover)  or  predators  (i.e,  present/absent  or 
fish/macroinvertebrates) . 

This  study  indicated  that  avoidance  of  sites  with  fish 
predators,  breeding  shortly  after  temporary  ponds  fill,  and 
rapid  development,  appear  to  be  the  major  anti-predator 
mechanisms  used  by  anuran  species  that  breed  in  temporary 
ponds,  and  that  within-wetiand  habitat  selection  was  of 
minimal  importance.  A better  understanding  of  the  potential 
role  of  microhabitat  use  as  an  anti-predator  strategy  could 
be  gained  by  examining  more  closely  the  size  distributions 
of  predator  and  prey  (e.g.,  Werner  et  al.  1983),  as  well  as 
the  foraging  strategies  of  predators  compared  to  the 
positional  use  by  tadpoles. 
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Table  3-4.  Correlations  between  all  species,  Hyla  squirella, 
Sana  utricularia,  and  Pseudacris  ocularis  with  water  depth, 
temperature,  pH,  plant  cover,  total  plant  biomass,  and 
macroinvertebrate  density  and  biomass. 
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[.’inure  3-4.  Variation  in  plant  cover  (mean  +1  SE)  among 
sampling  dates  (top),  wetland  (middle),  and  wetland  zone 
(bottom) . Zone  means  with  different  letters  are 
significantly  different. 
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CHAPTER  4 

EFFECTS  OF  COVER  AND  PREDATOR  SIZE  ON  SURVIVAL  AND 
DEVELOPMENT  OF  RANA  UTRICULAR  I A TADPOLES 

Introduction 

Size-limited  predation  is  a particularly  important 
process  during  amphibian  development.  Predation  levels  in 
systems  with  size-limited  predators  are  largely  dependent  on 
the  relative  sizes  of  predator  and  prey  (Ebenman  and  Persson 
1988)  . Numerous  studies  have  shown  that  predation  rates  on 
tadpoles  are  a function  of  tadpole  body  size,  and  a majority 
of  mortality  due  to  predation  occurs  early  in  development 
le.g.,  Banks  and  Beebee  19B8;  Cronin  and  Travis  1986; 
Richards  and  Bull  1990,  Semlitsch  1990,  Semlitsch  and 
Gibbons  1988;  Travis  et  al.  1985a).  Thus,  it  has  been 
suggested  that  selection  for  rapid  growth  during  the  larval 
phase  may,  at  least  in  part,  be  a response  to  predation  by 
size-limited  predators  (Travis  1983,  Travis  et  al.  1985b) . 

For  anurans  breeding  in  temporary  aquatic  habitats, 
oviposition  shortly  after  the  breeding  site  fills  provides  a 
potentially  important  mechanism  for  decreasing  predation 
pressure.  For  example,  because  dragonflies,  a major  predator 
of  tadpoles  in  temporary  breeding  sites,  also  oviposit  after 
pond  filling  (Ward  1992),  tadpoles  that  grow  rapidly  often 
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are  able  to  reach  a size  refuge  and  decrease  their 
cumulative  risk  of  predation.  Because  a majority  of 
predation  on  tadpoles  occurs  early  in  development,  the 
"growth  race"  between  dragonfly  naiads  and  tadpoles  can  have 
a particularly  large  influence  on  the  ability  of  predators 
to  regulate  tadpole  populations  {Caidwel!  et  al.  1980; 

Tejado  1993;  Travis  et  al.  1985a). 

Not  all  species  that  use  temporary  wetlands  employ  this 
reproductive  strategy,  however.  The  southern  leopard  frog 
(Rana  utricularia)  breeds  year  round;  however,  the  majority 
of  breeding  is  done  during  the  winter-spring  breeding 
season.  In  addition,  this  species  has  very  broad  breeding 
site  associations.  This  can  result  in  larvae  being  exposed 
to  an  assemblage  of  aquatic  insect  predators  that  are 
already  established,  and  therefore  large  in  size  relative  to 
newly  hatched  Ladpoles.  Thus,  early  predation  pressure  by 
large  predators  may  be  significant.  Under  such 
circumstances,  physical  features  of  the  environment,  such  as 
complex  habitat  structure,  may  provide  an  important 
mechanism  for  lowering  predation  rates. 

Habitat  structural  complexity  has  been  shown  to  be  an 
important  factor  mediating  predator-prey  interactions 
(reviewed  in  Heck  and  Crowder  1991).  Increased  habitat 
complexity  can  reduce  predation  rates  by  providing  cover  or 

Rozas  and  Odum  1988) , or  by  decreasing  foraging  success  of 
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predators  because  of  decreased  maneuverability  or  visual 
range  (e.g.,  Crowder  and  Cooper  1982;  Savino  and  Stein  1982; 
Werner  et  al.  1983).  R.  utricularia  tadpoles  can  grow  to  a 
large  size  relative  to  most  aquatic  insect  predators,  thus 
mechanisms  that  decrease  early  predation  pressure  may  have  a 
particularly  important  effect  on  overall  survival  rates. 

In  this  study  I exposed  tadpoles  of  the  southern 
leopard  frog  to  two  levels  of  cover  and  two  size  classes  of 
the  same  predator.  By  doing  so  I could  determine  whether 
increased  habitat  complexity  (i.e.,  cover)  provides  a 
mechanism  for  decreasing  predation  and  further  whether  such 
protection  is  dependent  on  predator  size.  1 predicted  that 
predation  rates  would  be  lower  at  high  cover  levels  and  that 
large  predators  would  be  more  effective  than  small  predators 
regardless  of  cover  level.  In  addition,  because  the  thinning 
effects  of  predators  can  release  prey  from  competition 
(Wilbur  1987,  1988),  I predicted  that  tadpoles  developing  in 
treatments  with  higher  predation  rates  (i.e.,  low  cover  and 
large  predator)  would  have  enhanced  growth. 

Methods 

Experimental  Design 

I examined  development  and  survival  of  tadpoles  in  a 
factorial  experiment  in  which  I manipulated  cover  level 
(high  versus  low)  and  predation  (large,  small,  or  no 
predator)  in  a complete  randomized  block  design  replicated 


four  times.  I established  a rectangular  array  of  24,  1.14  m 
diameter  plastic  wading  pools  with  tight  fitting  fiberglass 
screen  tops  at  the  MacArthur  Agro-Ecology  Research  Center 
(MAERC1  in  Lake  Placid,  Highlands  County,  Florida.  On  8 
March  1995,  I filled  pools  with  95  1 of  well  water  (depth  = 
12  cm)  and  added  1000  ml  of  well-raised  algae  and  zooplankton 
collected  from  several  wetlands.  I added  an  additional  500 
ml  of  algae  and  vegetation  (see  below)  on  16  March  1995. 
Water  overflow  was  prevented  by  a drainage  pipe.  I checked 
pools  weekly  to  determine  if  water  levels  were  maintained  at 
the  proper  level  by  rainfall.  When  rainfall  was  not  adequate 
to  compensate  for  evaporation,  1 added  well  water  (3  times 
during  the  experiment  I . 

I provided  two  levels  of  covet,  2000g  or  500g  (wet 
weight),  of  the  aquatic  plant  Hvdrochloa  carol inensis.  H. 
carolinensis  is  a thin-leafed  grass  that  provides  complex 
structure  throughout  the  water  column.  It  is  abundant  in 
wetlands  at  MAERC.  Both  treatments  provided  cover  throughout 
the  water  column;  however,  the  higher  cover  level  provided 

provided  relatively  sparse  cover.  Predator  treatments  were 
large  (mean  *_  1 SE:  4.18  1 0.26  cm;  n - 6)  or  small  (2,35  1 
0.11  cm;  n = 6)  'I'ramea  Carolina  (Odonata:  Libellulidae) 
naiads,  or  no  predator. 

On  23  March  1995,  I haphazardly  added  30  Gosner  (1960) 
stage  25  (0.04  • 0.01  g n = 15)  tadpoles  to  each 
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added  predators  the  following  day.  When  metamorphosis 
commenced,  I checked  pools  twice  a day  for  tadpoles  with 
emerged  forelimbs  (stage  42;  Gosner  1960) . I held  tadpoles 
individually  in  500-ml  plastic  jars  until  tail  resorption 
(stage  46;  Gosner  1960),  and  then  recorded  wet  mass  (to 
0.001  g) . 

During  the  second  week  of  May  1995,  the  area 
experienced  an  unseasonal  heat  wave  with  daytime  highs 
exceeding  32°C  for  several  days.  At  the  end  of  the  heat  wave 
air  temperatures  exceeded  38"C  and  water  temperatures 
reached  40“C.  This  event  caused  a massive  die-off  of 
tadpoles  on  15  May  1995.  At  this  stage  of  the  experiment  71". 
of  tadpoles  surviving  to  that  point  already  had 
metamorphosed  successfully.  I collected  the  dead  tadpoles 
and  sacrificed  any  tadpoles  that  remained  alive  and  obtained 
developmental  stage.  Thus,  15  May  1995  was  considered  the 
end  of  the  experiment.  Because  of  the  die-off,  response 
variables  were  the  mass  of  tadpoles  that  reached  Gosner 
(1960)  stage  46  (stage  42  by  15  May  1995),  age  of  metamorphs 
(calculated  as  the  start  of  the  experiment  to  stage  42),  and 
percent  survival.  Because  tadpoles  were  large  compared  to 
prudaLors,  I assumed  that  any  Ladpoles  in  Lite  pools  on  15 
May  1995,  whether  dead  or  alive,  would  have  survived 
predation.  Thus,  survival  was  based  on  the  number  of 
tadpoles  that  reached  stage  46  plus  pre-metamorphic  tadpoles 
that  remained  at  the  time  of  Lhe  die-off. 


Statistical  Analysis 


Treatment  effects  were  analyzed  using  two-way  fixed 
effect  analysis  of  variance  (ANOVAI . Mean  values  per  pool 
were  analyzed  because  measurements  of  individuals  within 
pools  are  not  independent.  One  replicate  from  the  no 
predator  x low  cover  treatment  was  dropped  from  the 
experiment  because  the  pool  developed  a leak  causing  rapid 
water  loss  resulting  in  tadpole  mortality.  Percentage  data 

log  transformed  (log(x  + 1)  for  mass  values)  prior  to 
analysis.  I performed  orthogonal  contrasts  to  test  for 
predictions  regarding  effect  of  predators  (high  predation  by 
large  predators)  and  cover  (higher  predation  at  lower 
cover),  as  well  as  prediction  regarding  enhanced  growth. 

When  interactions  were  significant  (p<0.05),  I conducted 
separate  comparisons  within  treatment  levels. 


Results 


Survival 

Predator  treatment  had  a significant  effect  on  total 
survival,  accounting  for  71. 6»  of  variation  in  survival. 


(Figure  4-1;  Table  4-1).  Orthogonal  contrasts  indicted  that 
survival  of  tadpoles  was  lower  in  treatments  with  large 
predators  compared  to  treatments  with  no  or  small  predators. 
This  was  true  at  high  (F=34.66,  pcO.OOl)  and  low  (F=11.80, 


interaction  between  predator 


SO 

and  cover  was  due  to  lower  survival  within  the  low  cover 
treatment  compared  to  the  high  cover  treatment  for  tadpoles 
exposed  to  large  predators.  However,  the  interaction 
explained  only  8.3S  of  variation.  Small  predators  were 
largely  ineffective.  Survival  of  tadpoles  in  pools  without 
predators  did  not  differ  from  that  of  tadpoles  in  pools  with 
small  predators,  regardless  of  cover  (contrasts:  low  cover 
p.0.11,  p=0.74;  high  cover  F=0.07,  p=0.80).  The  main  effect 
of  cover  was  not  significant.  Overall,  tadpole  survival  was 
lowest  under  the  large  predator  ::  low  cover  treatment. 

Mass  at  Metamorphosis 

Mass  at  metamorphosis  was  affected  by  the  main  effects 
of  predator  and  cover,  and  by  the  interaction  of  cover  and 
predator  (Figure  4-1;  Table  4-1).  Within  the  large  predator 
treatment,  tadpole  mass  was  higher  within  the  low  cover 
treatment  compared  to  the  high  cover  treatment  (F=24.42, 
p=0.002) . Effects  of  small  (F-0.02,  p=0.88)  and  no  predator 
( F=1 . 16,  p=0 . 30)  treatments  were  not  significantly 
influenced  by  cover.  Effects  of  predator  treatment  on 
tadpole  mass  at  metamorphosis  were  not  significantly 
different  within  the  high  cover  treatment  (P=2,748,  p=0.12); 
however,  in  low  cover  treatments,  tadpoles  within  the  large 
predator  treatment  were  significantly  heavier  than  those  in 
the  small  (F=37.51,  p<0.001)  and  no  (F*24.5,  p<0.001) 


predator  treatments.  Tadpole  mass  did  not  differ  between  the 
small  predator  and  no  predator  treatments  within  the  low 
cover  treatment  (F=0.dl , p=0.53). 


Age  at  Metamorphosis 

The  main  effect  of  predator  and  the  interaction  of 
predator  and  cover  had  a significant  effect  on  age  at 
metamorphosis;  however,  the  main  effect  of  cover  was  not 
significant  (Figure  4-1;  Table  4-1).  At  high  cover  levels, 
age  at  metamorphosis  was  not  significantly  affected  by 
predator  treatment  (F-0.01,  p=0.92) . However,  within  the  low 
cover  treatment,  tadpoles  exposed  to  large  predators 
metamorphosed  earlier  (F=11.34,  p<0.001)  compared  to  the 
small  and  no  predator  treatments,  which  did  not  differ 
(F=0.001,  p-0.97). 


Discussion 

Predator  sire  had  a much  larger  effect  on  tadpole 
performance  than  did  level  of  habitat  complexity.  High  cover 
level  did  reduce  predation  levels  compared  to  low  cover  when 
large  predators  were  present,  suggesting  that  habitat 
complexity  does  provide  increased  protection  from  predators. 
The  most  parsimonious  interpretation  for  lack  of  effects  by 
small  predators  is  that  tadpoles  were  able  to  grow  large 
enough  early  in  the  experiment  to  reach  a sire  refuge  and 
escape  predation.  Thus,  early  growth  by  R.  utricularia 
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tadpoles  developing  in  temporary  environments  can  greatly 
reduce  predation  pressure  from  aquatic  insects  that  are 
oviposited  around  the  same  time. 

Similarly,  growth  probably  decreased  overall  predation 
rates  by  large  predators;  however,  tadpoles  were  within  the 
gape  limits  o£  large  predators  for  a longer  time,  allowing 
increased  predation  by  odonate  naiads.  Because  tadpole 
survival  was  higher  at  high  cover  levels  when  large 
predators  were  present,  this  suggests  that  increased  cover 
decreased  the  foraging  efficiency  of  T.  Carolina  naiads.  A 
possible  alternative  hypothesis  is  that  tadpoles  raised 
under  higher  cover  grew  faster  and  therefore  decreased 
predator  success  by  surpassing  the  gape  limits  of  predators. 
This  is  unlikely,  however,  because  growth  of  tadpoles  on  the 
other  treatments  differed  little  relative  to  cover  level. 
Thus,  the  most  likely  explanation  for  decreased  predation  at 
high  cover  levels  is  interference  in  predator  foraging. 

Decreased  foraging  efficiency  with  increasing  habitat 
complexity  has  been  found  in  other  studies  of  aquatic  insect 
predators,  including  species  that  are  active  foragers  (e.g.. 
Anas  iunius)  and  sit  and  wait  predators  (e.g.,  Belostoma 
sp.)  (Chapter  6).  Heck  and  Crowder  (1991)  predicted  that 
less  mobile  predators  would  actually  have  increased  foraging 
efficiency  in  complex  habitats  due  to  increased  perch  sites 
for  foraging.  Banks  and  Beebee  (1988)  demonstrated  that 
increasing  complexity  from  no  vegetation  to  some  vegetation 


increases  foraging  success  of  odonate  naiads  on  Bufo 
calami ta  tadpoles.  Their  study  provides  some  support  for 
Heck  and  Crowder's  prediction.  However,  a majority  of 
research  to  date  suggests  that,  like  their  mobile  piscine 
counterparts,  aquatic  insects  often  have  decreased  foraging 
efficiency  in  structurally  comple::  habitats. 

R.  utricularia  tadpoles  reduce  activity  in  the  presence 
of  Ana:i  junius  naiads  (Chapter  5)  . Reduced  activity  has  been 
identified  as  a potentially  important  anti-predator 
mechanism  in  many  amphibian  species,  particularly  among 
species  that  breed  in  permanent  aquatic  habitats  where  fish 
are  the  dominant  predators  (Werner  and  McPeek  1994;  Werner 
and  Anholt  1993) . Research  suggests  that  the  anti-predator 
response  is  general,  rather  than  predator  specific  (Stauffer 
and  Semlitsch  1993) . Thus,  reduced  activity  may  have 
afforded  R.  ulriculatia  some  protection  from  predation.  The 
results  of  this  study,  however,  suggest  that  size 
differences  between  predator  and  tadpole  may  be  more 
important  than  prey  activity,  at  least  when  aquaLic  insects 
are  the  predators,  for  example,  increased  cover  also 
decreased  predation  levels  on  Hyla  squirella  tadpoles,  a 
species  LliaL  has  high  aeUvity  levels  even  In  the  presence 
of  predators  (Chapter  6) . 

The  thinning  effects  of  predators  resulted  in  larger 
tadpole  size  and  more  rapid  development  (earlier  age  at 
metamorphosis) . These  results  agree  with  other  studies  of 


predation  on  tadpoles  that  have  found  enhanced  growth  of 
surviving  tadpoles  (Wilbur  1987,  1988) . Under  these 
circumstances,  predation  can  provide  an  overall  positive 
effect  on  population  size,  resulting  in  a mutual istic 
relationship  between  predators  and  tadpoles  that  escape 
predation.  Particularly  in  environments  where  habitat 
desiccation  is  a significant  mortality  threat,  the  thinning 
effects  of  predators  can  lead  to  increased  survival  among 
individuals  that  escape  predation  (Wilbur  1987).  For 
example,  in  cattle  tank  experiments  where  habitat  drying  was 
manipulated,  Wilbur  (1987)  demonstrated  that  R . utricularia 
tadpoles  subject  to  predation  by  newts  actually  had  higher 
survival  compared  to  tadpoles  in  tanks  without  predators. 
When  predators  were  absent,  competitive  effects  led  to 
decreased  growth  and  prolonged  development,  resulting  in 
death  due  to  desiccation. 

In  this  study,  the  enhanced  growth  of  tadpoles  in  the 
large  predator  treatments  suggests  that  tadpoles  in 
treatments  where  predation  was  low  experienced  intraspecific 
competition.  Even  so,  this  had  little  effect  on 
susceptibility  to  predation,  which  was  determined  more 
strongly  by  piudatur  sire  than  any  competitive  etlecls  among 
tadpoles.  This  suggests  that  competitive  effects  were 
probably  not  significant  until  after  tadpoles  had  already 
grown  beyond  the  gape-limits  of  T.  carol ina  naiads. 
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The  importance  of  relative  size  differences  between 
tadpoles  and  aquatic  insect  predators  has  been  demonstrated 
many  times  (Brodie  and  Formanowicz  1983;  Caldwell  et  at. 
1980;  Crump  1984;  Formanowicz  19B6;  Richards  and  Bull  1990; 
Travis  et  al.  1985a).  Although  some  studies  have  shown  that 
larger  tadpoles  may  be  selected  over  small  tadpoles  ii  both 
size  classes  are  within  the  gape-limits  of  predators  (Crump 
and  Vaira  1991;  Tejedo  1993),  predation  on  tadpoles  that  are 
capable  of  growing  beyond  the  gape-limits  of  predators  is 
concentrated  on  smaller,  early-stage  tadpoles.  This  study 
provides  additional  evidence  of  the  importance  of  relative 
size  differences  between  predators  and  tadpoles.  Further,  it 
demonstrates  that  habitat  structural  complexity  can  play  an 
important  role  in  mediating  predator-prey  interactions,  even 
when  tadpoles  start  out  at  a size  disadvantage  relative  to 


predators. 
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Table  4-1.  Summary  of  ANOVA  for  responses  of  Rana 
utricularia  tadpoles  to  cover  and  predator  treatments.  CD  is 
the  coefficient  of  determination,  which  is  the  percentage  of 
variation  in  the  response  variable  attributable  to  the 
treatment  effect. 


Source  of  v 
Predator 

Predator  x Cover 


(1.03834  0.68  0.578 
1.19807  31.92  <0.001 
0.03493  1.86  0.194 
0.13963  3.72  0.051 
0.26271 


Source  of  variation 
Predator 


42.2 

14.3 
17.9 


Age  (days) 

Source  of  variation 

Block 
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0. 00009 
0.00088 
0.00007 
0.00107 
0.00153 
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Figure  4-1 . Mean  (i  1 SE)  responses  o£  Sana  utricularia 
tadpoles  to  cover  and  predator  treatments:  percent  survival 
(top),  wet  mass  at  metamorphosis  (middle),  and  age  at 
metamorphosis  (bottom) . Squares  indicate  high  cover,  circles 
indicate  low  cover. 


CHAPTER  5 

EFFECTS  OF  FOOD  AVAILABILITY  AND  RISK  OF  PREDATION  ON 
BEHAVIOR  AND  GROWTH  OF  RANA  UTRICULARIA  TADPOLES 


predators  are  broader  than  those  associated  with  direct 
predator-induced  mortality.  Studies  documenting  keystone 
predator  effects  (Paine  1969)  and  trophic  cascades 
(Carpenter  et  al.  1987;  Power  1990)  demonstrate  the  large 
effects  predators  can  have  on  assemblage  structure  even  when 

Direct,  but  non-lethal,  effects  (sensu  Strauss  19911 
can  have  a large  impact  on  prey  populations  as  well.  Such 
effects  are  often  manifested  through  a shift  in  the  behavior 
of  prey  in  response  to  the  presence  of  a predator. 

shift  in  microhabitat  use  or  a change  in  timing  or  amount  of 
activity  (reviewed  in  Lima  and  Dill  1990).  Facultative 
reductions  in  activity  in  response  to  the  presence  of  a 

including  aquatic  insects  (McPeek  1990;  Peckarsky  et  al. 
1993),  crustaceans  (Stein  and  Magnuson  1976),  fishes  (Fraser- 
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and  Gilliam  1987),  and  anuran  larvae  (Lawler  1989;  Semlitsch 
and  Reyer  1992;  Skelly  1999;  Skelly  and  Werner  1990;  Werner 
1991).  Reduction  in  activity  may  be  a widespread  response 
to  predators  because  prey  activity  is  a large  component  of 
the  process  through  which  many  predators  detect  prey  (Lima 
and  Dill  1990;  Werner  and  Anholt  1993).  Thus,  there  should 
be  strong  selection  pressure  against  activity  in  the  face  of 
predation. 

Reduction  in  activity,  however,  involves  a trade-off 
between  decreasing  risk  of  predation  and  obtaining  resources 
because  both  risk  (predation)  and  reward  (resources)  are 
positively  related  to  activity  (Werner  1992b;  Werner  and 
Anholt  1993) . The  trade-off  between  predation  risk  and 
resource  acquisition  suggests  that  an  important  non-lethal 
effect  of  predators  on  prey  is  a decrease  in  growth  (Skelly 
and  Werner  1990).  Both  intraspeeific  and  interspecific 
differences  in  activity  have  been  related  to  competitive 
ability  in  tadpoles:  active  tadpoles  are  often  competitively 
superior  (Lawler  1989;  Morin  1983;  Werner  1991,  1992a; 
Woodward  1983) . However,  active  tadpoles  are  also  more 
susceptible  to  predation  (Anholt  and  Werner  1995;  Azevedo- 
Ramos  et  al.  1992;  Skelly  1994).  This  inverse  relationship 
between  susceptibility  to  predation  and  competitive 
superiority  associated  with  activity  level  can  have  a large 
influence  on  assemblage  structure  (Morin  1983,  1986) . 


The  strength  of  non-lethal  predator  effects,  in  turn, 
may  be  largely  influenced  by  background  resource  levels. 
Theoretical  models  predict  that  activity  should  be  inversely 
related  to  mortality  risk  (Abrams  1993;  McNamara  and  Houston 
1987,  1994;  Werner  and  Anholt  1993).  These  models  also 
predict  that  foraging  activity  should  be  inversely  related 
to  resource  level  when  mortality  risk  is  dependent  on 
foraging  activity  (Abrams  1993;  McNamara  and  Houston  1987, 
1994;  Werner  and  Anholt  1993).  For  example,  Anholt  and 
Werner  (1995)  found  that  activity  rates  of  Rana  catesbeiana 
tadpoles  exposed  to  predators  were  inversely  related  to  food 
level,  and  this  resulted  in  lowered  mortality  for  tadpoles 
on  a high  food  level.  However,  when  mortality  risk  is 
independent  of  activity,  foraging  activity  should  increase 
with  increasing  resource  levels  (McNamara  and  Houston  1994) . 
Therefore,  both  predation  risk  and  resource  availability 
should  be  examined  to  understand  the  consequences  of 
behavioral  adjustments  on  individual  growth. 

The  purpose  of  this  study  was  to  examine  the 
behavioral  responses  and  growth  performance  of  tadpoles  of 
the  southern  leopard  frog  (Rana  utricularia)  raised  on 
different  resource  levels  both  in  the  presence  and  absence 
of  a non-lethal  predator.  By  altering  both  the  perceived 
risk  of  predation  and  resource  levels,  I was  able  to  examine 
the  relative  importance  of  these  factors  and  how  they 
interacted  during  the  development  of  anuran  larvae. 


Predictions 


Growth  and  Development 

The  Gilliam-Werner  model  (Gilliam  1982;  Werner  1986) 
predicts  that  fitness  in  pre-reproductive  individuals  is 
maximized  by  minimizing  the  ratio  of  mortality  (u)  to  growth 
rate  (gl . Based  on  this  model,  factors  that  increase 
mortality  in  the  larval  environment  should  result  in 
decreased  size  at  metamorphosis  (see  hypothetical  u/g  curves 
in  Werner  1986) . Increases  in  resource  level  increases 
growth  rate  and  therefore  results  in  the  opposite  effect; 
Increased  size  at  metamorphosis.  Reduction  in  activity  in 
response  to  predators  decreases  foraging  rates  and  therefore 
should  result  in  smaller  size  at  metamorphosis. 

The  Gilliam-Werner  model  makes  no  prediction  about  the 
length  of  the  metamorphic  period  because  reduced  size  could 
result  from  either  early  metamorphosis  at  smaller  size  or 
prolonged  larval  period  under  poorer  growth  opportunities 
resulting  in  smaller  size  (Wilbur  and  Collins  1973). 

If  tadpoles  are  responding  to  predators  they  should 
reduce  activity  in  their  presence  (Abrams  1993;  McNamara  and 
Houston  1987,  1994;  Werner  and  Anhoit  1983).  Further, 
tadpoles  under  a higher  food  regime  should  be  less  active 
than  those  under  a poor  regime  of  nutrition  in  the  presence 
of  a predator  (Werner  and  Anhoit  1993) . When  predators  are 


absent,  activity  rates  should  be  positively  related  to 
resource  level  (McNamara  and  Houston  1994). 

Based  on  these  models,  I predicted  that  tadpoles  would 
reduce  activity  in  response  to  the  presence  of  predators  and 
that  this  would  result  in  a smaller  size  at  metamorphosis. 
Further,  I predicted  that  activity  level  would  be  positively 
related  to  food  resource  when  predators  were  absent  but 
negatively  related  to  resource  level  when  predators  were 
present . 

Methods 

I examined  growth  and  behavior  of  tadpoles  in  a 
factorial  experiment  in  which  1 manipulated  resource  levels 
(low  versus  high  per  capita  food  level)  and  non-lethal 
exposure  to  larval  dragonflies  Ana;:  iunius  (Odonata: 
Aeshnidae) . On  19  February  1995  I collected  a single  R. 
utricularia  egg  mass  from  an  ephemeral  wetland  at  the 
MacArthur  Agro-Ecology  Research  Center  in  Lake  Placid, 
Florida.  The  egg  mass  was  held  in  a plastic  wading  pool 
filled  with  rain  water  and  allowed  to  hatch.  Late  instar  A. 
Iunius  naiads  were  collected  from  ephemeral  wetlands  at  the 
research  center  and  maintained  in  the  laboratory  on  a diet 
of  R.  utricularia  tadpoles. 

On  1 March  1995  I haphazardly  added  four  Gosner  (I960) 
stage  25  tadpoles  (mean  t 1 SE:  SVL  5.75  • 0.3  mm, 
preserved  wet  mass  = 0.031  * 0.003  g)  to  each  of  24  plastic 
containers  (34x29x14  cm).  Containers  were  filled  with  seven 
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liters  of  aged  well  water  and  arranged  on  laboratory  shelves 
in  six  blocks  of  four  containers  each.  Fiberglass  bags  (8x12 
cm)  were  suspended  along  one  side  of  the  container  to  hold 
predators.  Bags  prevented  predators  from  preying  on  tadpoles 
buL  allowed  tadpoles  to  chemically  detect  predators.  Bags 
in  predator  treatments  contained  one  late  instar  A.  iunius 
naiad,  whereas  empty  bags  were  used  as  controls  in  no- 
predator treatments.  Once  a week,  following  water  changes, 
each  predator  was  fed  a single  R.  utricularia  tadpole. 
Predators  were  replaced  with  new  individuals  as  necessary 
throughout  the  experiment. 

Tadpoles  were  fed  a finely  ground  mixture  of  Purina 
rabbit  chow  and  pelleted  fish  food  (3:1  by  mass).  Low  food 
treatments  received  per  capita  rations  (T.55  body  mass)  that 
were  limiting,  whereas  high  food  treatments  received  per 
capita  rations  (15*  body  mass)  that  promoted  growth  and 
development  (Alford  and  Harris  1988;  Werner  1992a).  The 
quantity  of  food  added  to  each  container  was  determined  by 
the  following  formula:  per  capita  food  ration  ::  mean  mass  of 
the  tadpoles  in  container  x number  of  tadpoles  surviving  in 
container  ::  number  of  days  to  next  feeding.  Within  each 
tood  treatment,  the  mean  tadpole  mass  per  container  was 
determined  once  a week  by  weighing  tadpoles  from  containers 
without  predators. 

I changed  container  water  every  week  at  the  beginning 
of  the  study,  twice  a week  by  week 


three,  and  every 
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day  by  week  five  to  avoid  water  fouling.  A light-dark  cycle 
of  14:10  was  provided  by  a series  of  tlorescent  lights  and 
natural  lighting.  Water  temperature  varied  between  22-24"C 
during  the  experiment. 

Behavioral  Observations 

1 used  scan  sampling  to  measured  activity  rates  and 
location  of  tadpoles  within  containers  during  the  first  four 
weeks  of  the  experiment  (n=18  observations) . I made 
observations  between  0900-1400  h on  non-feeding  days  and 
terminated  observations  at  4 weeks  because  older  tadpoles 
appeared  lo  react  to  my  presence.  1 recorded  the  number  of 
tadpoles  moving  (swimming  or  feeding)  and  the  number  of 
tadpoles  on  the  side  of  the  container  opposite  the  predator 
bag  and  converted  counts  to  percentages  for  each  bin 
averaged  across  the  18  observation  periods. 

Response  Variables  and  Statistics 

Metamorph  mass  (wet  weight)  was  measured  when  at  least 
one  forelimb  had  emerged  (stage  42,  Gosner  1960) . The  larval 
period  was  defined  as  the  start  of  the  experiment  to 
forelimb  emergence.  Container  means  were  used  as  response 
variables  for  all  statistical  analyses.  Hass  and  larval 
period  data  were  log-transformed  and  activity  and 
distribution  data  were  angularly  transformed  to  homogenize 
variances  and  normalize  data  prior  to  analysis.  I used 
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factorial  ANOVA  testing  for  Lhe  effects  of  laboratory  block, 
predator  presence,  resource  level,  and  the  interaction  of 
predator  presence  and  resource  level  on  response  variables 
(Snedecor  and  Cochran  1980) . I performed  orthogonal 
contrasts  to  test  my  predictions  regarding  effects  on  mass 
and  activity.  Additional  contrasts  determined  whether  the 
presence  of  a predator  affected  the  length  of  the  larval 
period. 


Besul ts 

Effects  on  Survival 

Survival  to  forelimb  emergence  was  981.  One  tadpole 
jumped  out  of  its  container  and  one  died  of  unknown  causes. 

Effects  on  Growth  and  Development 

Food  treatment  had  a significant  effect  on  tadpole  mass 
at  metamorphosis,  accounting  for  491  of  the  variation  in 
tadpole  mass  (Table  5-1).  Tadpoles  on  the  high  food 
treatment  were  25. VI  larger  than  tadpoles  on  the  low  food 
treatment.  The  main  effect  of  predator  was  not  significant; 
however,  there  was  a significant  interaction  between  food 
and  predator  treatment  (Figure  5-1).  Tadpoles  on  the  low 
food  treatment  were  14.31  smaller  when  predators  were 
absent,  whereas  tadpoles  on  the  high  food  treatment  were 
8.21  larger  in  the  absence  of  A.  Junius.  Orthogonal 
contrasts  indicated  that  differences  in  mass  within  food 
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treatment  were  significant  for  tadpoles  on  low  food  (F=4.88, 
p=0.043)  but  not  for  tadpoles  on  high  food  (F»1.659, 


Both  food  and  predator  treatments  affected  age  at 
metamorphosis  (Table  5-1).  Food  had  the  largest  effect, 
explaining  11.1%  of  the  variation  in  larval  period,  whereas 
predator  treatment  accounted  for  only  8.51  of  the  variation 
(Figure  5-1).  Lowered  food  levels  resulted  in  a 12.31 
increase  in  larval  period,  whereas  non-lethal  exposure  to  a 
predator  prolonged  the  larval  period  by  4.31  (about  2 days). 
The  interaction  between  food  and  predator  was  marginally 
significant  (Table  5-1).  Orthogonal  contrasts  showed  that 
tadpoles  on  low  food  reached  metamorphosis  significantly 
earlier  when  predators  were  absent  (F-14.15/  p=0,002); 
however,  age  at  metamorphosis  was  not  affected  by  predator 
treatment  for  tadpoles  on  high  food  (F=2.32,  p=0.148). 


Effects  on  Tadpole  Activity  and  Distribution 

The  presence  of  A.  iunius  had  a significant  effect  on 
the  spatial  distribution  of  tadpoles,  accounting  tor  45.2". 
of  the  variation  in  distribution  (Table  5-2).  Overall,  69.8% 


ol  the  tadpoles 
the  predator  whe 
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sn  A.  iunius  was  present  compared  to  54.81 
vere  absent  (Figure  5-2) . Neither  food  nor 
of  food  and  predator  was  significant. 
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The  effects  on  activity  wete  more  complex.  Predator 
presence  had  a significant  effect;  tadpoles  in  predator 
containers  reduced  activity  by  60.81  (Figure  5-2).  Doth  food 
and  the  interaction  of  predator  and  food  treatments  were 
marginally  significant  (Table  5-2) . When  predators  were 
present,  activity  levels  were  similar  on  both  food 
treatments.  Activity  levels  were  higher  when  predators  were 
absent  for  both  the  high  food  (F«29.57,  p=0.0007)  and  low 
food  (F=5.94,  p=0.028)  treatments;  however,  tadpoles  on  the 
high  food  treatment  had  a higher  level  of  activity  compared 
to  tadpoles  on  the  low  food  treatment  when  predators  were 
absent  (Figure  5-2). 


Discussion 

Behavioral  Responses 

Tadpoles  altered  both  spatial  distribution  and  activity 
levels  in  response  to  the  presence  of  predators.  Because 
activity  level  is  positively  correlated  with  mortality  risk 
ISkelly  1995;  Werner  and  Anholt  1993),  reduced  activity  can 
be  an  important  mechanism  for  decreasing  risk  of  predation. 
Reductions  in  activity  in  response  to  predators  have  been 
Lound  in  several  other  aiiuraii  species  including  Pseudaeris 
crucifer  (Skelly  1995),  Bufo  terrestris  (Skelly  and  Werner 
1990),  Rana  clamitans  and  Rana  catesbeiana  (Werner  1991, 
1992b),  Rana  lessonae  and  Rana  esculenta  (Horat  and 
Semlitsch  1994) . Decreases  in  tadpole  activity  in  response 
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to  predators  appears  to  be  a general,  rather  than  a 
predator-specific  response  (Lawler  1989;  Semlitsch  and 
Gavasso  1992).  Further,  interspecific  differences  in 
response  to  predators  appear  to  be  strongly  related  to 
differences  in  breeding  habitat  association.  Species  that 
behaviorally  respond  to  predators,  such  as  those  listed 
above,  typically  develop  in  sites  that  contain  large 
predator  populations  (i.e.,  permanent  sites  or  sites  with 
long  hydroperiods) . In  contrast,  species  using  the  most 
ephemeral  habitats,  which  often  lack  predators,  generally  do 
not  possess  the  same  behavioral  mechanisms  for  avoiding 
predation  (Babbitt,  unpubl.  data;  Woodward  1 183). 

The  costs  and  benefits  of  behavioral  activity  patterns 
in  tadpoles,  and  therefore  the  expression  of  such  patterns, 
appear  to  be  strongly  related  to  differences  in  the  major 
source  of  mortality  at  permanent  (i.e.,  predation)  versus 
ephemeral  sites  (i.e.,  habitat  desiccation) (Skelly  1995). 

In  most  situations,  behavioral  adjustments  that  are  flexible 
and  context-specific  should  be  strongly  selected  for  because 
individuals  possessing  such  flexibility  should  be  better 
able  to  balance  costs  and  benefits  associated  with  specific 

that  are  too  short  in  duration  for  larvae  to  complete 
development,  species  with  tadpoles  that  reduce  mortality 
risk  through  reductions  in  activity  could  potentially  breed 
in  a wide  variety  of  sites  containing  different  suites  of 


99 

predators.  R.  utricuXarla  provides  a good  example,  because 
it  breeds  in  a wide  range  of  habitats,  including  permanent 
sites  containing  fish  predators  and  temporary  sites  where 
macroinvertebrates  are  the  main  predators  (pers.  obsl. 

R.  utricularia  tadpoles  also  altered  spatial 
distribution  in  response  to  predators.  In  simple 
environments,  such  as  the  one  used  in  this  study,  spatial 
shifts  are  limited  to  near  versus  away  from  the  predator 
(e.g.,  Skelly  and  Werner  1990).  Other  studies  have 
demonstrated  shifts  to  refuge  areas  and  benthic  micohabitats 
(Holomuzki  1986b;  Lawler  1989;  Semlitsch  and  Gavasso  1992). 
The  potential  growth  costs  of  shifts  in  spatial  distribution 
are  clear  for  species  that  occupy  refuge  areas  that  lack 
food,  or  for  species  that  shift  locations  to  areas  with 
lower  food  availability.  However,  for  individuals  that  are 
developing  in  areas  that  offer  partial  refuges  from 
predators,  such  as  areas  with  increased  plant  cover, 
movement  away  from  predators  may  not  always  Involve  a trade- 
off between  mortality  and  growth  (Diehl  1992) . Increased 
cover  may  provide  increased  food  resources,  as  well  as 
increased  protection  from  predators  (Mclvor  and  Odum  1988) . 
Overall,  the  largest  costs  associated  with  predator- induced 
habitat  shifts  are  probably  linked  to  differences  in  the 
spatial  distribution  of  food  resources,  competitors,  and 
predators.  Thus,  trade-offs  between  mortality  risk  and  food 
acquisition  related  to  habitat  shifts  may  or  may  not  be 


biologically  significant. 


Growth  and  Development 


The  presence  of  predators  resulted  in  longer  larval 
periods  (i.e.,  delayed  metamorphosis)  on  both  food 
treatments.  However,  the  effects  on  mass  at  metamorphosis 
were  dependent  on  resource  level.  Tadpoles  raised  under  the 
low  food  treatment  metamorphosed  at  larger  sizes  when 
predators  were  present,  whereas  metamorphs  on  the  high  food 
treatments  were  larger  when  predators  were  absent.  Although 
tadpoles  had  similar  activity  levels  when  predators  were 
present,  tadpoles  on  the  high  food  treatment  were  more 
active  when  predators  were  absent  compared  to  tadpoles  on 
the  low  food  treatment.  This  suggests  that  tadpoles 
subjected  to  predators  on  the  high  food  treatment  incurred  a 
higher  cost  in  terms  of  growth  due  to  a greater  decrease  in 
activity. 

Thus,  although  high  food  availability  resulted  in 
larger  size  compared  to  low  food  availability  regardless  of 
predator  treatment,  the  cost  of  reducing  activity  in 
response  to  predators  was  higher  on  high  food.  A possible 
explanation  for  lack  of  a similar  cost  for  tadpoles  on  the 
low  food  treatment  is  that  because  less  food  was  available, 
even  at  lower  activity  levels  tadpoles  were  able  to  garner 
all  resources  in  the  container  within  each  feeding  period, 
whereas  tadpoles  on  the  high  food  predator  present 
treatment  were  not.  In  a similar  study,  Skelly  and  Werner 
(1990)  found  no  interaction  of  predator  and  food  treatment 


for  either  age  or  mass  at  metamorphosis.  In  their  study* 
decreased  activity  associated  with  predator  presence  did  not 
affect  larval  period,  but  resulted  in  reduced  size  at 
metamorphosis.  In  field  enclosures,  Skelly  (1992)  found  that 
Hyla  versicolor  tadpoles  responding  to  predators  had  both 
decreased  growth  and  developmental  rates.  Thus,  although 
costs  of  predator-induced  behavioral  adjustments  appear  to 
be  common,  how  these  costs  are  reflected  in  growth  versus 
development  appear  to  vary  among  tana  or  under  varying 
experimental  conditions  (Werner  and  Anholt  19911. 

Several  studies  have  found  that  increases  in  food  level 


are  associated  with  decreased  activity  (reviewed  in  Werner 
and  Anholt  1993) . Such  patterns  have  been  found  in  Rana 
svlvatica  (Werner  1992a),  Rana  catesbeiana  (Anholt  and 
Werner  1995),  and  Pseudacris  triseriata  (Skelly  1995).  In 
contrast,  Bufo  terrestris  (Skelly  and  Werner  1990)  and 
Pseudacris  crucifer  (Skelly  1995)  showed  little  responses  to 
food  levels.  Kohler  and  McPeek  (1989)  found  higher  activity 
levels  at  lower  food  levels  in  a Baetis  mayfly  and  a 
Glossosma  caddisfly.  In  this  study,  R.  utricularia  had 


true  for  activity  responses  to  predators,  differences  among 
anuran  species  in  their  responses  to  resource  level  may  be 
related  to  breeding  habitat  associations,  and  the  relative 


of  reducing  activity  (Skelly  1995). 
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predict  that  species  breeding  in  permanent  ponds  would  be 
more  likely  to  decease  activity  at  high  resource  levels 
because  the  costs  of  decreased  growth  are  comparatively 
small.  In  contrast,  for  species  developing  In  more  ephemeral 
sites,  reduced  foraging  activity  could  greatly  increase  the 
possibility  of  death  due  to  habitat  drying.  For  these 
species,  we  would  predict  either  small  adjustments  or 
increased  foraging  activity  at  higher  resource  levels 
(Skelly  1995)  . 

Ecological  Consequences  of  Responses 

The  duration  of  the  larval  period  has  important 
consequences  for  species  that  use  temporary  habitats  as 
breeding  sites.  Because  the  threat  of  mortality  due  to 
habitat  drying  is  significant,  factors  that  prolong  the 
larval  period  can  have  negative  consequences  (Skelly  1995) . 
In  this  study,  tadpoles  responded  to  predators  by  prolonging 
the  larval  period.  R.  utricularla  breeds  in  a variety  of 
habitats,  ranging  from  temporary  wetlands  to  permanent 
ponds.  Depending  on  specific  breeding  habitat,  the 
prolonged  larval  period  may  or  may  not  increase  risk  of 
inorLaliLy  due  Lo  liabiLat  deslccaliuii.  Ueducliuns  in  acLivily 
can  have  positive  effects  by  decreasing  risk  of  mortality 
regardless  of  habitat  association  because  for  most  of  Lhe 
broad  suite  of  predators  that  R.  utricularla  may  face  in 


different  breeding 


whether  vertebrate 
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invertebrate,  prey  movement  is  an  important  component  of 
predator  detection  (Lima  and  Dili  1990) . Therefore,  the 
costs  of  reduced  activity  in  terms  of  prolonged  larval 
period  should  be  weighed  against  the  benefit  of  decreased 
predation  risk.  Because  R.  utricularia  breeds  in  temporary 
sites,  it  would  be  instructive  to  lest  the  response  of 
tadpoles  to  the  presence  of  predators  under  varying  drying 
regimes  to  determine  whether  R.  utricularia  can  balance  the 
costs  and  benefits  of  foraging  versus  predator  detection 
with  the  costs  of  habitat  drying. 

Reductions  in  size  at  metamorphosis  can  have  potential 
costs  regardless  of  breeding  habitat  association.  Size  at 
metamorphosis  can  have  an  important  influence  on  adult 
survival  and  reproductive  fitness  (Berven  1990;  Goater  1994; 
Semlitsch  et  al.  1988;  Smith  1987).  In  this  study,  tadpoles 
on  the  higher  plane  of  resources  that  responded  to  predators 
incurred  relatively  higher  cost  in  terms  of  reduced  size  at 
metamorphosis,  although  they  had  absolutely  larger  masses 
compared  to  tadpoles  on  the  lower  food  treatment.  Food 
treatment  had  a greater  effect  on  growth  than  did  predator 
treatment.  Increased  size  at  metamorphosis  as  a response  to 

in  Wilbur  1980) . That  food  availability  had  such  a large 
effect  on  growth  is  not  surprising;  however,  the 
differential  costs  of  responding  to  predators  on  different 
food  treatments  is.  These  results  demonstrate  the  need  for 


resource  level  and  threat  of 


assessing  responses  to  both 
predation  in  order  to  fully  understand  the  consequences  of 
cost-benefit  trade-offs  mediated  by  changes  in  behavior. 

Interspecific  differences  in  behavioral  responses  to 
predators,  and  the  consequences  these  differences  have  on 
growth  and  survival,  are  important  considerations  for 
explaining  assemblage  structure  at  different  breeding  sites 
or  habitat  patches  (Werner  and  Anholt  1996) . Across 
breeding  sites,  interspecific  differences  in  behavioral 
responses  to  predators  largely  may  represent  evolutionary 
responses  to  different  costs  associated  with  breeding 
habitats.  For  example,  the  relative  mortality  costs  due  to 
predation  versus  habitat  drying,  may  largely  explain 
interspecific  differences  in  behavioral  responses  to 
predators  (Shelly  1995;  Woodward  1983).  Within  breeding 
sites,  the  effects  of  interspecific  differences  in 
behavioral  responses  on  growth,  coupled  with  patchy 
distribution  patterns  of  predators  and  tadpoles,  could 
influence  local  differences  (i.e.,  among  patches)  in 
assemblage  structure.  Results  of  empirical  studies  (Anholt 
and  Werner  1995;  Werner  and  Anholt  1996)  as  well  as 
theoretical  models  (MacNamara  and  Houston  1987)  suggest  that 
population  regulation  must  be  viewed  as  a multi-factor 
process.  Behavioral  responses  made  by  individuals  to  changes 


acquisition, 


therefore  population-level  responses. 


Whether  behavioral  adjustments  have  direct  or  indirect 
effects  within  food  web  dynamics,  our  understanding  of 
processes  will  be  enhanced  by  considering  behavior  as  i 
important  factor  in  future  experiments- 
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CHAPTER  6 

EFFECTS  OF  COVER  AND  PREDATOR  IDENTITY  ON  PREDATION  OF 
HYLA  SQU1RELLA  TADPOLES 

Introduction 

Habitat  structural  complexity  can  strongly  affect 
biotic  interactions.  In  particular,  the  intensity  of 
predator-prey  interactions  can  be  mediated  by  habitat 
complexity  because  the  effectiveness  of  predators  often 
decreases  in  structurally  complex  habitats  (Babbitt  and 
Jordan  1996;  Diehl  1988;  Werner  et  al.  1983).  Typically, 
increased  habitat  complexity  reduces  predation  rates  by 
providing  refuges  for  prey  or  by  decreasing  predator 
efficiency  (Crowder  and  Cooper  1982;  Savlno  and  Stein  1982; 
Werner  et  al.  1983). 

The  role  of  habitat  structure  in  mediating  predator- 
tadpole  interactions  is  potentially  important  because 
predation  is  a significant  source  of  mortality  of  anuran 
tadpoles  (Calef  1973;  Beyer  el  al.  1979;  Smith  1983). 

Species  breeding  in  permanent  ponds  display  a wide  array  of 
anti-predator  defenses  for  avoidinq  predation  by  fish 
(Formanowicz  and  Brodie  1982;  Hats  et  al.  1988;  Lawler 
1989).  In  general,  such  defenses  are  lacking  in  species  that 
breed  in  ephemeral  ponds  (Kats  et  al.  1980).  In  addition, 
activity  rates  of  tadpoles  in  temporary  ponds  often  are 
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higher  compared  to  species  in  permanent  ponds  (Woodward 
1983) . Because  many  predators  of  tadpoles  are  gape-limited, 
rapid  growth  to  a size-refuge  is  an  important  mechanism  for 
escaping  predation  (Crump  1984 ; Richards  and  Bull  1990; 
Travis  et  al.  1985a;  but  see  Crump  and  Vaira  1991).  For 
species  breeding  in  ephemeral  environments,  rapid  growth  is 
also  important  for  decreasing  risk  of  desiccation  (Newman 
1989;  Wilbur  1987).  However,  trade-offs  exist  between  rapid 
growth  and  predator  avoidance  because  actively  foraging 
tadpoles  are  more  likely  to  be  detected  by  predators  (Skelly 
1994;  Werner  and  Anholt  1993;  Woodward  1983).  Although 
ephemeral  ponds  lack  fish  predators,  these  sites  often 
contain  aquatic  insect  predators  that  also  are  significant 
mortality  agents  (Smith  1983;  Wilbur  and  Fauth  1990) . 
Therefore,  any  mechanism  that  can  reduce  predator  efficiency 
without  decreasing  foraging  activity  could  enhance  tadpole 
survival,  particularly  in  ephemeral  environments.  Because 
structurally  complex  habitats  can  reduce  predator  foraging 
efficiency,  complex  environments  may  provide  tadpoles  with 
partial  protection  from  predators,  even  when  tadpoles  are 
actively  foraging. 

In  this  study  1 examine  how  habltaL  structure  allecls 
predation  rates  on  squirrel  treefrog  (Hyla  souirella) 
tadpoles.  H.  squirella  is  a common  anuran  in  Florida,  and  is 
often  numerically  dominant  in  the  ephemeral  sites  in  which 
it  breeds  (Babbitt,  unpubl.  data).  H.  squire! la  tadpoles  are 


active  foragers  even  in  the  presence  of  insect  predators 
(Babbitt,  unpubl.  data),  and  therefore  they  are  a good  model 
to  test  whether  increased  habitat  structure  provides 
protection  for  actively  foraging  tadpoles. 

Methods 

I conducted  the  experiment  at  the  MacArthur  Agro- 
Ecology  Research  Center  in  Lake  Placid,  Highlands  County, 
Florida.  1 established  a rectangular  array  of  16,  1.14  in 
diameter  plastic  wading  pools  with  tight-fitting  screen 
tops.  On  3 September  1994  I filled  pools  with  95  1 of  well 
water  (depth  = 12  cm) . To  allow  normal  feeding  behavior,  I 
added  2000  ml  of  pond  water  from  a wading  pool  used  for 
culturing  algae.  Pools  received  either  1000  g (high  cover) 
or  500  g (low  cover)  of  the  wetland  plant  Hvdrochloa 
carolinensis.  This  plant  is  common  in  many  wetlands  where  H. 
squire! la  breeds,  and  provides  complex  structure  throughout 
the  water  column. 

1 collected  one  gravid  female  and  several  males  and 
placed  them  together  in  an  aquarium  with  aged  well  water  and 
sparse  vegetation.  The  female  layed  a clutch  of  eggs  on  1 
September  1994,  and  on  15  September  1994  1 haphazardly 
assigned  30  Gosner  (1960)  stage  35-37  tadpoles  (mean  * 1 SD; 
SVL  12.2+0.47  mm,  weight  0.301+0.04  g;  n-25)  to  each  pool. 
Prior  to  the  experiment,  1 kept  tadpoles  together  in  a 
plastic  wading  pool  and  provided  rabbit  chow  ad  libitum.  I 
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used  late-stage  tadpoles  because  Hie  predators  I used  are 
extremely  effective  predators  of  early  stage  tadpoles 
(Babbitt,  unpubl . data),  and  therefore  using  early  stage 
tadpoles  would  not  provide  a good  test  of  habitat  effects. 

I used  larval  common  green  darners  (Ana::  junius, 
Odonata:  Aeshnidae)  (total  length  29.3^0.14  mm;  n=8)  or 
giant  water  bugs  (Lethocerus  americanus.  Hemiptera; 
Belostomatidae)  (total  length  17. MO. 09  nun;  n=8)  as 
predators.  Both  species  are  common  inhabitants  of  the 
ephemeral  wetlands  in  which  H.  souirella  breed.  1 held 
predators  individually  without  food  for  48  hrs  to 
standardize  hunger  within  each  species. 

1 began  the  experiment  at  1900  hrs  on  16  September  1994 
when  1 added  a single  predator  to  each  pool.  I terminated 
the  experiment  at  1900  hrs  on  22  September  and  censused 
remaining  tadpoles  to  determine  survival.  Diurnal  high 
temperature  was  31"C  during  the  experiment. 

The  experiment  was  a 2X2  factorial  design  arranged  in  a 
randomized  complete  block  and  replicated  four  times.  Two 
species  of  predator  (A.  junuis  or  t.  americanus)  were 
crossed  with  two  plant  densities  (1000  g or  500  g) . In 
addition  to  predaLor  treatments,  1 also  ran  two  control 
pools  at  each  plant  density  without  predators  to  determine 
whether  plant  treatment  alone  affected  survival.  I used 
fixed  effect  analysis  of  variance  (ANOVA)  to  determine 
whether  predator  identity,  plant  density,  or  their 


interaction  affected  survival  of  H.  squirei  la  tadpoles.  The 
number  of  tadpoles  surviving  was  log  transformed  to  reduce 
skewness  (Snedecor  and  Cochran  1980) . Controls  were  not 
included  in  the  above  analysis  because  all  tadpoles  survived 
in  the  control  pools,  and  all  mortality  in  treatment  pools 
could  be  attributed  to  predation. 

Results 

Tadpoles  had  higher  survival  in  the  high  cover 
treatment  compared  to  low  cover  (F:ii,-26.15;  p=0.003). 
Predator  identity  did  not  significantly  affect  survival 
(F,  , =1.63;  p=0.26),  and  the  interaction  between  predator 
species  and  plant  density  was  not  significant  { F.t  =0.64; 
p=0.44) . The  number  of  tadpoles  surviving  under  the  higher 
density  treatment  was  nearly  double  that  of  the  low  density 
treatment  under  both  predator  treatments  (Figure  6-1). 

Discussion 

Compared  to  research  on  many  other  aquatic  organisms, 
only  limited  empirical  research  has  focused  on  effects  of 
habitat  complexity  in  mediating  relations  between  predators 
and  anuran  larvae.  Sredl  and  Collins  (1992)  found  that 
differences  in  habitat  structure  did  not  affect  survival  of 
mountain  treefrog  (Hvla  enimia)  tadpoles  subject  to 
predation  by  tiger  salamanders  (Ambystoma  tiqrinum) . 
Similarly,  Figiel  and  Semlitsch  (1991)  examined  predation  by 
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crayfish  (P.  acutus)  on  injured  gray  treefrog  (H. 
chrysocelis)  tadpoles  and  found  no  effects  of  habitat 
complexity.  In  contrast,  increased  plant  density  resulted 
in  decreased  predation  by  aquatic  insects  on  southern  toad 
IBufo  terrestris)  tadpoles  (Babbitt  and  Jordan  1996). 
Differences  among  these  studies  could  be  related  to  the 
relative  degree  of  cover,  or  to  differences  in  predator  or 
prey  behavior. 

Reduction  in  foraging  success  for  the  aquatic  insects 
used  in  this  study  probably  resulted  from  both  reduced 
visual  range  and  interference  in  stalking  behavior.  For 
example,  visual  detection  is  important  for  prey  capture  by 
odonate  naiads  (Pritchard  19651,  and  increased  vegetation 
cover  can  interfere  with  the  predatory  efficiency  of  these 
active  foragers  (Diehl  1992;  Folsom  and  Collins  1984) . 
However,  foraging  efficiency  of  L.  americanus,  a sit  and 
wait  predator,  also  was  reduced  suggesting  that  increased 
habitat  complexity  can  mediate  predator-prey  interactions 
between  tadpoles  and  a broad  suite  of  aquatic  insect 
predators.  For  example,  in  a similar  study,  predation  on 
B,  terrestris  tadpoles  was  lower  at  higher  cover  levels  for 
both  odonate  naiads  and  hemipLeians  (Belostoma  f luminea) 
(Babbitt  and  Jordan  1996)  . 

H.  sguirella  tadpoles  develop  rapidly  and  are  active 
foragers,  even  in  the  presence  of  aquatic  insect  predators 
(Babbitt,  unpubl.  data).  Rapid  growth  and  development. 


necessary  for  attaining  metamorphlc  size  prior  to  pond 
drying,  may  explain  why  many  species  developing  in  ephemeral 
sites  do  not  respond  to  predators.  Although  active  tadpoles 
are  generally  more  susceptible  to  predators  (Skelly  1994; 
Werner  and  Anholt  1993;  Woodward  1983),  species  developing 
in  ephemeral  ponds  presumably  would  be  more  restricted  in 
reduced-activity  responses  If  the  threat  of  habitat  drying 
poses  a more  significant  mortality  threat.  Rather,  foraging 
in  structurally  complex  environments  may  provide  increased 
protection  from  predators,  even  for  species  that  do  not 
facultatively  decrease  activity  in  response  to  predators. 
Research  that  examines  the  relative  role  of  habitat 
structure  in  mediating  predation  by  different  predator 
groups  (e.g.,  fish  versus  aquatic  insects)  on  species  with 
different  activity  rates  should  provide  valuable  insight 
into  the  relative  costs  and  benefits  of  different  behavioral 
strategies  and  the  degree  to  which  habitat  features  can 
influence  their  effectiveness. 


Low  High 

Cover  Level 


Figure  6-1.  Number  of  Hyla  squirella  tadpoles 
predation  under  different  cover  levels.  Bars  i 


CHAPTER  7 

SUMMARY  AND  CONCLUSIONS 


The  Introduction  of  water  from  fish-containing  ditches 
into  otherwise  isolated  temporary  wetlands  produced  large 
changes  in  species  composition  of  anurans.  Wetlands  Lhat 
experienced  the  hydrologic  disturbance  were  avoided  by 
species  that  breed  in  fish-free  sites;  however,  use  by 
species  that  do  breed  in  sites  with  fish  increased.  Species 
that  avoided  impacted  sites  bred  in  adjacent  sites  that  were 
unaffected  by  water  from  ditches.  These  results  suggest  that 
temporary  wetlands  on  this  site  provide  dynamic  habitats 
that  offer  varying  breeding  opportunities  and  conditions  for 
larval  developmental  that  are  highly  dependent  on 
meteorological  conditions.  The  effect  of  the  hydioiogic 
disturbance  is  a system  that  is  both  spatially  and 
temporally  dynamic  resulting  in  variable  assemblage 
structure  and  composition  at  affected  sites.  Because  MAEP.C 
contains  numerous  wetlands  and  ditches  that  vary  in 
hydrology  from  ephemeral  to  permanent,  most  species  can  find 
suitable  breeding  habitat  except  under  extreme  drought 
conditions.  Ditching  exacerbates  the  drought  by  increasing 
drainage,  lowering  the  water  table,  and  prolonging  dry 
conditions. 
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Historical  changes  Lo  the  landscape  have  probably 
increased  local  anuran  species  richness  and  altered  tadpole 
assemblage  structure  at  the  wetlands  at  MAERC.  The 
surrounding  upland  matrix  can  have  a large  influence  on 
wetland  use  as  breeding  siLes.  At  least  one  species,  llvla 
femoral is,  appeared  to  be  locally  abundant  but  limited  in 
distribution  to  forested  hammocks.  Research  on  the 
relationship  between  upland  habitat  patch  characteristics 
and  wetland  use  would  increase  our  understanding  of  breeding 
use  patterns,  as  well  as  our  ability  to  provide  guidance  for 
managing  anuran  populations  in  altered  systems. 

By  examining  spatial  and  temporal  patterns  of  tadpoles 
and  aquatic  insect  predators,  I was  able  to  determine  that 
microhabitat  use  within  wetlands  does  not  appear  to  be  an 
important  mechanism  for  avoiding  predation  within  the 
temporary  wetlands  that  I studied.  The  invertebrate 
predators  that  are  the  dominant  predators  at  these  sites 
were  relatively  evenly  distributed  throughout  these 
wetlands;  therefore,  differential  use  of  wetland  areas  by 
tadpoles  presumably  would  not  have  provided  increased 
protection.  Examinations  of  microhabitat  use  by  Hyla 
squirel  la,  a specie's  Hint  bleeds  sltoilly  after  wetland 
filling,  and  Rana  utricularia  which  is  more  plastic  in 
breeding  site  use,  revealed  only  weak  microhabitat 
selection.  Early  breeding,  rapid  development,  and  avoidance 
of  sites  containing  fish  predators  appear  to  be  the  major 
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squirella. 


contrast. 


utricularia  may  rely  more  heavily  on  early  growth  (as 
opposed  to  rapid  development)  to  attain  a sire  that  is 
beyond  the  gape  limits  of  predators,  and  behavioral 
adjustments  in  activity  to  avoid  detection  by  predators. 
Such  differences  in  life  history  strategies  influence  the 
range  of  breeding  sites  that  each  species  can  use 
successfully. 

Some  studies  have  demonstrated  spatial  differences  in 
tadpole  distribution  within  wetlands.  Most  of  these  studies 
have  been  conducted  in  wetlands  with  zonal  differences  that 
vary  greatly  in  cover,  such  as  wetlands  with  open-water 
centers.  The  differences  in  cover  among  zones  in  the 
wetlands  I studied  were  reiatively  small,  and  may  not  have 
provided  enough  difference  in  protection  or  food  (or  other 
factors)  to  generate  differences  in  microhabitat  use. 
Tadpoles  may  exhibit  plasticity  in  microhabital  selection 
relative  to  the  degree  of  difference  in  microhabitat 
features  such  as  cover.  Interspecific  differences  in 
tadpole  plasticity  relative  to  cover  could  lead  to 
differences  in  assemblage  structure  under  varying  habitat 
conditions  (i.e.,  dill-:  • -iti-Lactoi 

experiments  that  examine  microhabitat  effects  are  needed  to 
determine  the  relative  importance  of  factors,  such  as  cover, 
in  structuring  tadpole  assemblage. 
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Although  tadpoles  did  not  differentially  occupy 
different  areas  of  wetlands,  experimental  research 
identified  the  importance  of  increased  cover  in  mediating 
predator-prey  interactions.  This  was  true  for  tadpoles  that 
behaviorally  adjust  activity  to  decrease  detection  by 
predators  (R.  utrlcularial  and  those  that  do  not  (H. 
squirella) . Further,  cover  decreased  foraging  rates  of 
aquatic  insects  that  are  sit-and-wait  predators  that  cling 
to  vegetation  (Lethocerus  americanus) , those  that  cling  to 
vegetation  but  are  more  active  (Ana::  iunius) , and  those  that 
are  benthic  Cl'raroea  Carolina) . Thus,  although  tadpole 
microhabitat  selection  relative  to  cover  did  not  appear  to 
be  important,  occupation  of  areas  with  high  cover  may 
decrease  predation  rates.  Higher  densities  of  aquatic 
insects  in  areas  with  high  cover  could  negate  the  benefits 
of  cover;  however,  in  this  study  predators  were  distributed 
evenly.  One  aspect  not  examined  in  this  study  was 
microhabitat  use  relative  to  the  water  column.  Differences 
in  positional  use  within  the  water  column  (i.e.,  mid-water 
versus  benthic)  could  be  an  important  anti-predator 
mechanism  for  some  species.  For  example,  a large  number  of 
benthic  predators  may  pose  liLtle  threat  to  tadpoles 
foraging  in  the  water  column. 

Finally,  tadpoles  of  some  species,  including  R. 
utricularia,  exhibit  shifts  in  spatial  location  and  activity 
levels  in  response  to  the  presence  of  predators.  Reduction 
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in  activity,  however,  involves  a trade-off  between 
decreasing  risk  of  predation  and  obtaining  resources  because 
both  risk  (predation)  and  reward  (resources)  are  positively 
related  to  activity.  However,  the  relative  costs  of 
decreasing  activity  in  response  to  predators  are  dependent 
on  background  resource  level.  In  R.  utricularia.  the 
presence  of  predators  resulted  in  longer  larval  periods 
(i.e.,  deiayed  metamorphosis)  regardless  of  food  treatment. 
However,  the  effects  on  mass  at  metamorphosis  were  dependent 
on  resource  level.  Tadpoles  raised  under  the  low  food 
treatment  metamorphosed  at  larger  sizes  when  predators  were 
present,  whereas  metamorphs  on  the  high  food  treatments  were 
larger  when  predators  were  absent. 

These  results  support  the  idea  that  population 
regulation  is  a multi-factor  process.  Behavioral  responses 
made  by  individuals  to  changes  in  the  environment  affect 
risk  of  predation,  resources  acquisition,  and  therefore 
population  level  responses.  Because  the  effects  of 
predation  can  be  direct  or  indirect,  our  understanding  of 
these  processes  will  be  enhanced  by  considering  behavior  as 
an  important  factor  in  future  experiments.  Because  long 
hydroperiod  temporary  wetlands  oiler  suitable  breeding 
habitat  to  a wide  variety  of  anurans,  the  potential  for 
behavioral  differences  among  the  species  developing  within 
these  sites  to  affect  assemblage  dynamics  may  be  great. 
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